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^ ■ ABSTRACT 
^' 

Qh' High-dispersion spectra centered on the Li 6708 A hne have been obtained for 70 

O ■ potential members of the intermediate-age open cluster NGC 3680, with an emphasis 

> , on stars in the turnoff region of the cluster CMD. A measurable Li abundance has been 
c/3 , 

, derived for 53 stars, 39 of which have radial velocities and proper motions consistent with 

I cluster membership. After being transferred to common temperature and abundance 

I scales, previous Li estimates have been combined to generate a sample of 49 members, 

' 40 of which bracket the cluster Li-dip. Spectroscopic elemental analysis of 8 giants and 

^ ■ 5 turnoff stars produces [Fe/H] = -0.17 ± 0.07 (sd) and -0.07 ± 0.02 (sd), respectively. 

^Sj ■ We also report measurements of Ca, Si and Ni which are consistent with scaled-solar 

QQ . ratios within the errors. Adopting [Fe/H] = -0.08 (§3.6), isochrone comparisons 

^ , lead to an age of 1.75 ± 0.1 Gyr and an apparent modulus of (m — M) = 10.30 ± 

I 0.15 for the cluster, placing the center of the Li-dip at 1.35 ± 0.03 Mq. Among the 

^ I giants, five of nine cluster members are now known to have measurable Li with A(Li) 
near 1.0. A combined sample of dwarfs in the Hyades and Praesepe is used to delineate 

^ • the Li-dip profile at 0.7 Gyr and [Fe/H] = +0.15, establishing its center at 1.42 it 0.02 

— ' Mq and noting the possible existence of secondary dip on its red boundary. When 



evolved to the typical age of the clusters NGC 752 (age = 1.45 Gyr, (m — M) = 8.4), 
IC 4651 (age = 1.5 Gyr, (m-M) = 10.4), and NGC 3680, the Hyades/Praesepe Li-dip 
profile reproduces the observed morphology of the combined Li-dip within the CMD's 
of the intermediate-age clusters while implying a metallicity dependence for the central 
mass of the Li-dip given by M/Mq = 1.38 ±0.04 + 0.4 ±0.2 [Fe/H]. The implications 
of the similarity of the Li-dichotomy among giants in NGC 752 and IC 4651 and the 
disagreement with the pattern among NGC 3680 giants are discussed. 



■^Visiting Astronomer, Cerro Tololo Inter-American Observatory. CTIO is operated by AURA, Inc., under contract 
to the National Science Foundation. 
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Subject headings: open clusters and associations:individual (NGC 752, NGC 3680, IC 
4651, Praesepe, Hyades) - stars: abundances - techniques: spectroscopic 



1. INTRODUCTION 



Open clusters are studied for a variety of reaso ns, ranging from their intended use as probes of 



the st ructural and chemical evolution of the Galaxy (jFriel et al. l2nn2l : lYong et al. 11200,4 



20081) to delineation of the f iner details of stellar structure and evolution (jMeibom et al. II2002I : 



Sestito et al. 



VandenBerg &: Stetson Il2004l l . While lithium is an element with clear cosmological origins and ram- 
ifications, studies of lithium in open clusters predominantly emphasize its role within stellar evolu- 
tion, though insights gained from these inv estigations could have an i ndirect impact upon our under- 



stand ing of the cosmological parameters (jSestito &: Randich 



2005 



Korn et al. 



2006 



Delivannis 



2000|). Lithium, boron and beryllium are relatively fragile nuclei. Attriti on of lithium from the pre - 
sumed initial value of A(Li) = 3.30 implied by meteoritic abundances (j Anders &: Grevesse Ill989l ) 
and very young star clusters, is expected in the internal layers of stars or protostars where the 
temperatures exceed ~ 2.5 x 10® K, above which Li burns via proton capture to ^He. Over a 
star's lifetime, the integrated effect of any structural features that put surface layers in contact 
with the deeper layers at or above this temperature will be expressed in the depletion of surface 
Li. The surface Li abui idance can also be enhanced if radiative acceleration pushes it i nto the 
surface convection zone (jRicher &: Michaud 1 119931 : iDelivannis. Steinhauer &: Jeffries 1 120021') or if a 



deep e ning convection zone d redges up previously diffused Li ([Delivannis. Demarque &: Kawaler 
1990; ISills Dellivanislboool ). 



StandarcQ stellar evolution theory (SSET) predicts the attrition of Li as a function of stellar 
mass, age, and composition. Low mass pre-main sequence models evolving down the Hayashi 
track develop high enough temperatures at the base of the surface convection zone (SCZ) so 
that energetic protons destroy Li there, thereby depleting the surface Li abundance. Lower stel- 
lar mass implies a deeper and longer-lived SCZ, which implies more Li depletion; there is lit- 
tle Li depletion for models reaching the ZAMS as A, F, and G dwarfs , while Li depletion is 
more severe for K and cooler spectral types (e.g. iPinsonneault I (|l997l )). Li depletion ceases 



as the SCZ recedes toward the surface and its base cools. Independently of the details of the 



MS for AFG dwarfs (e.g. ( 


Proffitt k Michaud 


1989: 


Delivannis. Demaraue &; Kawaler 


1990; 


Pinsonneault. Kawaler &: Demarque 


1990 : Swenson et al. 


I994I: Chabover et al. 


I995I1) whereas 



KM dwarfs can continue to deplete Li as they evolve beyond the ZAMS. The predicted Li-Te rela- 
tion is modified very slightly as models evolve to slightly higher Tg over time. Higher-metallicity 
models have deeper convection and deplete more Li. 



^By 'standard' we mean spherically symmetric models that ignore effects due to diffusion, rotation, mass loss and 
magnetic fields. 
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These predictions are qualitatively consistent with observations of main sequence stars in 



youn ger open clusters. Observations of stars in young clusters like the Pleiades (jSoderblom et al 



19931 ) do confirm the general ZAMS Li-Tg relation resulting from pre-main-sequence evolution. 
However, stars in real clusters exhibit major departures from the theoretical expectations. A sig- 
nificant dispersion in the Li abundance is found at any given stellar mass or surface temperature, 
implying that factors beyond the star's basic structural parameters of mass, chemical composition 



and age are relevant (jKing et al. II2OOC1 ). An additional clue is posed by the observation that short- 



period, tidally-locked binaries retain their high initial surface Li, an association that highlights 



rotation and/or stellar rotational hi story as additional factors (Deliv annis et al. 1 1 19941), in agree- 



199S 



Mallik et al. 



2003; 



Pasquini et al. II2004I ). 



Jasniewicz et al. 



ment with studies of both field stars (iRandich et al. Ill999l : IChen et al. , 2001.: .do Nascimento et al. 
2OO3I:' " ' 



20061 ) and clusters (jThorburn et al. 



1993 



Jones et al. 



Among older open clusters, depletion of Li in low mass stars as a function of stellar mass 
exceeds SSET predictions by amounts that grow increasingly disparate with age. At ages between 
0.1 to 0.2 Gyr, a feature strikingly at odds with SSET begins to appear in p lots of Li abundance 
versus stellar surface temperature or mass (jSteinhauer &: Delivannis II2004I ). By the age of the 
Hyades (0.65 Gyrs), stars with surface temperatures of 6600ib200 K, consistent with main sequence 
masses near 1.4 Mq, show abrupt depletions in Li of up t o 2 orders of magnitude, defining the Li 
or Boesgaard-Tripicco dip (IBoesgaard &: Tripicco Ill986al ). Such stars would not be expected to 
have convective envelopes; no obvious physical mechanism exists to put the surface layer in contact 
with the hotter, deeper layers. As emphasized by the many papers evaluating the properties 



relevant to this dist i nct feature (|Pasauini 



et al. 



2001 



Chen et al. 



2001 



Mallik et al. 2003; Pasquini et al. 2004; Steinhauer Sz Deliyannis 20041 ). mapping the depth of 



Boesgaard fc King II2OO2I : 



the dip and its range in temperature and mass as a function of metallicity and age is critical to 



any a ttempt to decipher the underlying physical mechanism that produces it (jCharbonnel &: Talon 



2OO5I ). Moreover, among clusters in the 3 to 5 Gyr range, the stars leaving the main sequence for the 
subgiant and giant branch emerge from the Li-dip with inherently lower Li that should be driven 
even lower by post-main sequence evolution. Discovery of any evolved stars with measurable Li in 
clusters of this age range should be a clear indicator of anomalous evolution, mass transfer, a.nd/or 
deep mixing of newly processed Li to the surface of stars on th e giant branch (jBalachandran Ill995l : 



Pasquini et al. 



2001 



2004 



Anthony- Twarog &: Twarog II2004I ). Finally, the common age of cluster 



dwarfs offers the opportunity of defining reliably the time evolution of the Li depletion among 
the cooler stars beyond the Li-dip if significant samples with sufficient s pectroscopic accuracy 



can be compiled over the a ge range of interest (jSestito Sz Randich I l2005l : iRandich et al. 
Prisinzano &: Randich 1120071 ). 



2006 



The focus of this investigation is the intermediate-age open cluster, NGC 3680. It was selected 
as a prelude to the analysis of a significantly more challenging but exceptional cluster, NGC 6253, 
for several reasons. As discussed in more detail in Sec. 3.3, some Li d ata for NGC 3680 stars from 
the turnoff through the giant branch do exist (jPasquini et al. II2OOII ) (PRP), though the sample 
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at and below the turnofF is by no means complete; we hoped to broaden the characterization of 
the Li-dip profile by expanding the sample with comparably high signal-to-noise spectra. At the 
same time, we wanted adequate overlap with previous work, not just so that conclusions might be 
strengthened by combined results for common stars in NGC 3680 but also to assist in validating our 
results for NGC 6253, based on spectra obtained with the same equipment and similarly reduced. 
Third, while photometric metallicity estimates are available for the turnoff stars and giants, as 
well as moderate-dispersion spectroscopic abundances for a few giants, metal ab undances based 



on high-dispersion spectroscopy have been obta ined for only a fe w cluster giants (IPasauini et al. 



20011 : ISmilianic et al. II2008I : ISantos et al. 1 120091 ') and two dwarfs (|Pace et al. II2008I : ISantos et al. 
2OO9I ') . one of which might not be a member. Obtaining sufficiently high S/N spectra for NGC 



3680's dwarf stars does not present an insuperable challenge, but the task is complicated by the 
severe paucity of single dwarfs in the cluster. 

Broader astrophysical circumstances make NGC 3680 an intriguing target for observation. It 
is one of a trio of well-studied clusters of similar age: NGC 752 is closer than NGC 3680 but of 
comparable metallicity, while IC 4651 is more richly populated, mor e distant and definitely more 



meta l -rich than eit h er cluster, as confirmed by both p hotometrv (lAnthonv-Twarog Twarog 



1987; 



Nissen 



1988; 



(jPasquini et al. I l2004l : 



Anthonv- Twarog Twarog 



Carretta et al. 



2004 



200c ; 



Meibom et al. 



Pace et al. 



2OO2I ) and spectroscopy 



20081 ) . However, the modest metallic- 



ity range of ~0.3 dex among the trio brack ets the range occupie d by the large majorit y of open 



clust e rs of any age withiri 1.5 kpc of the Sun ([Twarog et al. 1119971 ). Along with NGC 752 ([Twarog 



19831 : [Daniel et al. Ill994ll. the color-magnitude diagram (C MP) of NGC 3680 exhibits a bimodal 



main sequence ([Nissen Ill988l : lAnthony- Twarog et al. Ill989l ). i.e., a secondary sequence displaced 



to the red of the comparably populated ma in sequence, the byproduct of a rich binary population 



(Daniel et al. 



1994 



Nordstrom 



Nordstrom et al 



et al 



19971 ). Having collated photometry, proper motions and ra- 
([19971 ) derived a binary fraction of 57% among -^100 members of 



dial velocities, 
NGC 3680 with complete data. 

NGC 3680, like IC 4651 and NGC 752, appears to have a poor to modest population below 
IM^ ([Anthonv-Twarog et al. Ill99ll : iHawley et al. Ill999l ). Based on the extant stellar population, 
Nordstrom et al. ( 19971 ) estimate that the cluster mass might have been ten times greater and the 
stellar population 30 times more numerous before its apparent evaporation. 

The intermediate age of ~ 1.7 Gyr for the three clusters places their turnoff stars squarely 
in the region of the Li-dip, but young enough that stars higher in mass than the dip should 
still be present near the turnoff region. The turnoff morphology also implies a turnoff mass 
that simultaneous ly corresponds to the rn a ss for which convective core overshoot should p lay a 



conspicuous role (Mazzei Pieatto 



1988 



Daniel et al. lll994l : lKozhurina-Platais et al. 



Anthonv- Twarog et al. I Il989l : lAndersen et al. ][l990; 
I997I ). Confirmation of these characteristics has been 



based on exquisite information a.b out proper- motion and radial- velocity niembe rships and binarity 



([Kozhurina-Platais et al. Ill995l : iMermilliod et al. Ill995l : iNordstrom et al.lll996l ). Once the numer 



ous binaries are identified and removed from the CMD, the location of the single-star main sequence 
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most resembles stellar models incorporating significant convective overshoot (jNordstrom et al 



19971 ) 



The layout of the remainder of the paper is as follows: Sec. 2 details the spectroscopic obser- 
vations and their reduction to produce radial velocities used in conjunction with proper motions 
to identify cluster members. Sec. 3 presents the fundamental parameters that play a role in the 
analysis of the spectra, as well as the derivation and interpretation of the abundances. Sec. 4 
discusses the Li abundance pattern in NGC 3680, NGC 752, IC 4651, the Hyades, and Praesepe. 
The composite structure of the Li-dip and its possible evolution are revisited in Sec. 5 and our 
conclusions are summarized in Sec. 6. 



2. SPECTROSCOPIC DATA 

2.1. Observations and Reductions 

Spectra of 71 stars in NGC 3680 were obtained at the Cerro Tololo Inter- American Observatory 
Blanco 4m telescope using the Hydra multi-object spectrograph. Three exposures for each of 
two fiber-position configurations were obtained on 18 and 19 July, 2005. One fiber configuration 
emphasized stars above the turnoff for a total exposure of one hour; the second reached further down 
the main sequence with a total exposure of four hours. Hydra was used in echelle mode; operating 
at relatively low order (8, in our case), dedicated filters select and pass the desired spectral region 
(filter 6 for our purposes). Our spectra were characterized by a dispersion of 0.15 A per pixel in 
the range 6520 A to 6810 A with S/N per pixel that ranged from ~160 to 20. Measured widths for 
calibration lamp lines imply a spectral resolution of 15, 000 in the region near 6700 A. 



The numbering system of the proper-motion study bv iKozhurina-Platais et al. I (|l995l )(KGP) 
is one of the larger sets; it will be used in the discussions of our own data. All exposures for 
star KGP 1396 were severely compromised by placement of the spectra on a region of the chip 
affected by a filter defect, so no further discussion of results for this star based on our data will be 
included. Processing of the data and calibration frames was completed at Indiana University and 
the University of Kansas using standard programs within IRAF. For clarity, all quoted uncertainties 
refer to the standard deviation unless explicitly noted otherwise. 



2.2. Radial Velocities 

Radial- velocity and rotational-velocity estimates were derived for 70 stars in our program 
utilizing the fourier-transform, cross-correlation facility, FXCOR, in IRAF. Program stars are com- 
pared to stellar templates of similar temperature, in this case over the wavelength range from 6575 
A to 6790 A including Ha, as well as a narrower region in the vicinity of Ha alone. A comparison 
of the results from these two internal measures demonstrates that, if two stars with significant 
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discrepancies are excluded from the comparison (one is very cool, the other has very low S/N), 
the average difference in the radial- velocity determinations, in the sense (full spectrum - Ha), is 
—1.27 ± 2.9 km/sec. The radial velocity for KGP 802 that appears in Table 1 is based on the Hq 
region alone; radial velocities for all other stars are based upon the full spectrum. 

To refer the derived radial velocities to an absolute scale, two options exist. First, two radial- 
velocity standards were observed on each of the nights for the CTIO Hydra runs, HR 9014 and 
HR 9032. A comparison of our determinations with the values from the General Catalog of Radial 
Velocities (jWilson Ill953l ) suggests that our velocities are too large b y 1.1 ±2.1 km/sec. Seco nd, we 
can also compare our radial velocities for 18 st ars in common with iNordstrom et al.l (|l996l ) and 3 
stars in common with iMermilliod et al. I (|l995l ). With known or suspected spectroscopic binaries 
excluded from the comparison, the average velocity difference, in the sense (N 96 - this study) 
is —0.45 ± 1.66 km/sec, indicating that our radial velocities are consistent with INordstrom et al 



(119961 ) within the errors. A similar cqnipari son to three apparently single-star red giants in common 
with the survey of IMermilliod et al. I (|l995l ) produces a velocity difference, in the sense (M95 - this 
study), of —0.75 ± 0.53 km/sec. We elected not to correct our radial velocity values based on 
these fairly noisy comparisons; additional discussion in the next session will suggest that our radial 
velocity values are larger than previous results by 0.4 km/sec. 



2.3. Cluster Membership 



Radial- velocity studies by IMermilliod et al. I (jl995l ) and INordstrom et al.l (jl996l ) established 
membership credentials for 15 red giants and over 100 red g iant and dwarf stars, respectively, 
followed by an a nalysis of the cluster's evolutio nary history (INordstrom et al.l Il997l ) . A nearly 



pa rallel study bv iKozhurina-Platais et al. I (119971 ) followed similar lines. The proper-motion study 



by Kozhurina-Platais et al. ( 19951 ) yields membership probabilities for 2711 stars. By themselves, 
the proper-motion measures do not discriminate cleanly between cluster members and field stars 
for NGC 3680. This point is illustrated in Fig. 1, a histog ram of the number of star s as a function 
of proper-motion probability for the 391 stars described bv lKozhurina-Platais et al. I (|l995l ) as "the 
most probable cluster members." From this figure alone, it's difficult to find a sharp distinction 
between members and i ion-members, except to conclud e that star s with P(iJ,) < 1 to 20% ar e 
unlikely to be members. iKozhurina-Platais et al. I (|l995l ). as well as INordstrom et al.l (|l996l . Il997l ). 
combine information from radial velocities and proper motions to establish more comprehensive 
membership criteria. In the discussion of membersh ips for our candidate star s that follows, we 
will employ the proper-motion based probabilities of iKozhurina-Platais et al. I (|l995l ) and radial 
velocities determined from our spectra. 

Fig. 2 shows our radial velocities, along with proper-motion-based membership probabilities 
for our 70 program stars, indicating a cluster mean radial velocity near 2 km/sec and illustrating the 
eventual criteria for cluster membership. Whether or not binaries are included, the median radial 
velocity for stars with P(/i) > 20% is 1.6 km/sec. The average radial velocity for 24 stars with 
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proper-motion membership probability P > 20% that are not known or suspected spectroscopic 
binaries, is 1.64 it 1.22 km/sec. If giants are excluded, the mean becomes 1.72 it 1.22 km/sec. 
The average radial velocity for the four single-star, red giant members (KGP 1175, 1374, 1461 and 
1873) is 1.28 lb 0.56 km/sec. As expected from the earlier discussion, our ru e an ra dial velocity 
from single member dwarfs is 0.45 km/sec larger than that of lNordstrom et al.l (|l997l ). who find a 
mean of 1.27 km/sec from dwarf members. Our ave rage based on red giants alone is similarly 0.41 
km /sec larger than the value 0.87 km/sec found by iMermilliod et al. I (jl995l ) from giant members. 
Nordstrom et al.l ()1997l ) comment on the difference between the average radial velocity of the main 
sequence and giant stars, and conclude that it is consistent with an expectation based on different 
surface gravities for the two classes of stars; we also find a difference of ~ 0.4 km/sec between the 
dwarfs and giant s , but the significance of this difference appears to be only marginal in our case. 
Nordstrom et al.l (|l997l ) estimate the internal velocity dispersion from dwarf stars as 0.65 km/sec; 
using an estimate of our internal accuracy of ±0.9 km/sec leads to a velocity dispersion of 0.8 
km/sec. 

En tries in Table 1 include th e more common identifications: KGP, as used by SIM BAD, refer- 
ring to iKozhurina-Platais et al. I (|l995l ) ; the well-known Eggen numbers (jEggen I Il969l ) ; WEBDA 
numbers, which use Eggen numbers for the brighter stars and refer to other surveys for the fainter 
stars. Photome tric V magnitudes and B — V colors are taken from the collated and merged values 
as described by I Anthony- Twarog &: Twarog I (|2004| )(ATT). Estimates of rotational velocity and ra- 
dial velocity follow, each with error estimates from FXCOR in IRAF. The final column summarizes 
information about binarity. Comments about binary status are as specific as possible, includ i ng nq - 
tations "SB?" referring to the analysis of radial- velocity errors by iKozhurina-Platais et al. I (119971 ) 
that suggest binary status and, similarly, "w/sb?" indicating that our spectra suggest unusually 
wide and/or doubled lines in the instrumental spectra. Stars with proper- motion probabilities and 
radial velocities within the bounds illustrated in Fig. 2, between —1.8 and 5.1 km/sec, are consid- 
ered members; note that while using a modest radial-velocity range to constrain membership, we 
have extended the proper- motion probability limit to 10%. 



3. Cluster Properties - NGC 3680 

3.1. Age and Distance Estimates 

ATT collated and mapped Stromgren h — y colors from four other surveys to their own photom- 
et ry in NGC 3680, then mapped th e magnitudes and colors to the broad-band photometric system 



of iKozhurina-Platais et al. I (|l997l ). These are the photometric values we will use in the analyses 
of our spectroscopic data below. ATT presented independent determinations of the interdependent 
quantities of foreground reddening and overall metal abundance for NGC 3680, as well as extensive 
reviews of earlier determinations of both; the reader is referred to this discussion for details. In 
summary, NGC 3680 was found to have E(i? — V) = 0.058 ± 0.003 (sem), well within the range 
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bounded by the UBV-hased estimate of 0. 04 bvlEggen I (11969) and the high-end deterni ination of 
0.10 from DDO photometry of the giants (jClaria &: Lapasset Ill983l : iTwarog et al. Ill997l ). 



Photometric metahicity estimates from the intermediate-band mi and hk indices led to a 
value of [Fe/H] = -0.14 ifc 0.03 dex (sem), lower than the earlier estiin ates of -1-0.09 to -1-0.11 



of Nissen 



(lAnthonv-Twarog et al. Ill989l : iNordstrom et al.lll997l : iBruntt et al. 1 1 19991 ) tied to the photometry 



( 19881 ). The lower metallicity removed the apparent discrepancy between the photo- 
metric dwarf abundances and thos e of the giants defined by DDO photometry, moderate-dispersion 
spectroscopy (jTwarog et al. 1119971 ). and high dispersion spectroscopy (PRP). Mindful that spectro- 
scopic analyses are necessarily tied to stellar parameters obtained assuming an a priori metallicity 
estimate, we adopted an initial value of [Fe/H] = —0.14 based upon the photometric results of ATT 
but iteratively settled on a value of [Fe/H] = —0.08. The rationale behind this choice is discussed 
in Sec. 3.4 - 3.6. 

With reddening and abundance estimates in hand, a comparison of the homogeneous photom- 
etry to appropriate isochrones constrains the cluster age and distance, while providing measures of 
the individual stellar surface gravities. The CMD based on the composit e BV photometric values 



(ATT ) for 63 probable members, selected using the joint proper-motion (jKozhurina-Platais et al 



19951 ) and radial-velocity criteria discussed in Sec. 2.3 from the stars observed for this investigation 
and those of PRP, is shown in Fig. 3, with open symbols denoting the 22 stars for which some binary 
status is known or suspected. Triangles are stars for which only the Li data of PRP are available. 
While the photometric data set is not substantively different from that used by ATT, the isochrones 
and interpolation software for specific ages and abundances ado pted for the pre s ent discussion are 



those of Y2 ( http:// www.astro.yale.edu/demarque/yyiso.htmll) (jYi et al. II2003I : iDemarque et al. 



20041 ). rathe r than the updated P adova isochrones ( http : / / stev .oap d . inaf . it / cgi-bin / cmd [) . The lat 



ter models (jMarigo et al. 1 12008| ) are substantially the same as the iGirardi et al. I (|2002l ) models 



used in ATT and require similar offsets of -1-0.01 mag and -0.023 mag, respectively, to the My and 



B — V color scales to bring a one Mq star of age 4.6 Gyr (jBahcall et al. 



1995 



Houdek fc Gough 



20081 ) into alignment with our consistently adopted values of 4.84 and 0.65 for the solar My and 
{B — V) color. The comparable shifts for the set are —0.03 and -1-0.008, respectively. The 
Padova solar isochrones are defined using the Z=0.019 set as the solar composition, while Y^ define 
solar as (Y,Z) = (0.266, 0.0181). The switch in isochrone sets allowed a comparison with earlier 
work to get an indication of how sensitive the conclusions are to the specific models adopted, but 
was explicitly done because the Y^ isochrones generate masses for the stars in the Li-dip that are 
more internally consistent with age, a point we will return to in Sec. 4 and 5. For Z=0.0152 ([Fe/H] 
= —0.08), Y^ isochrones with ages of 1.7 Gyr and 1.8 Gyr are displayed in Fig. 3, adjusted for the 
adopted reddening of E(i3 — V) = 0.058 and an apparent distance modulus of (m — M) = 10.3. 
The age and apparent distance modulus are in good agreement with those discussed in ATT (1.85 
Gyr, 10.2); the revised Padova isochrones yield 1.85 Gyr and (m — M) = 10.2, while supplying a 
slightly better match to the shape of the turnoff below the hook. 



Note in Fig. 3 that the stars traditionally identified as the clump have the predicted luminosity 
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for He-core-burning stars in this phase but overlap the colors of the first-ascent giant branch. 
Moreover, the luminosity of the majority of the stars also is consistent with the location of the 
red giant bump, where first-ascent giants temporarily reverse their evolution up the giant branch 
when the H-burning shell crosses the discontinuity in composition created by the deep limit of the 
convective atmosphere as the star ascends the giant branch. 

The solid agreement between observation and theory in Fig. 3 boosts confidence in the derived 
stellar parameters needed for the spectroscopic analysis that can be deduced differentially from the 
CMD. As exemplified by the Padova comparison, however, even assuming a fixed reddening and 
metallicity, adoption of a different set of isochrones based upon alternate models could result in 
slightly different age, distance, surface gravity and mass estimates for the individual stars. From 
the detailed discussion in ATT, the spread in (m — M) found in the recent literature for NGC 
3680 is a byproduct of the combined variation in the derived [Fe/H], the adopted reddening, and 
th e choice of stellar i sochro nes, as exemplified by the valu e s of ( m — M) = 10.45 to 10.65 derived 
by iNordstrom et al.l (119971 ) and iKozhurina-Platais et al. I (jl997l ) adopting [Fe/H] much closer to 
Hyades metallicity for the cluster, as well as alternate sets of isochrones. 



3.2. Cluster Metallicity - Li 



SPLOT, part of the one-dimensional spectrum analysis suite in IRAF, was used to measure 
equivalent widths of the line blend incorporating the Li-doublet at 6707.8 A and the Fe I line 
at 6707.45 A; the EW measures are summarized in Table 2, along with the estimated S/N for 
each spectrum. S/N per pixel estimates were based on the modal statistics of the spectra in a 
"relatively" line-free region between 6615 A and 6625 A. Li measurements were not possible in 
some stars: either the line was too weak to be measured or the spectrum was distorted in a way 
which made measurement implausible. Of the eight stars in the latter category, three have S/N per 
pixel below 30 (KGP 802, 803 and 1610), while the others (KGP 979, 988, 1347, 1410A and 1506) 
are known or suspected spectroscopic binaries whose multiple components confused the Li-line 
region beyond reclamation. 

The largest set of Li EW measurements and abundance estimates for NGC 3680 is that of PRP, 
offering abundances for 26 stars, including 2 suspe cted non-members, ari d expanding on an earlier 
discussion of four approximately solar- type stars by lRandich et al. I (|2000l ). While the larger sample 
of PRP will supply the comparison of choice for our measurements, the discrepancy between the 
published equivalent widths from the same authors for the four stars based upon the same spectra 
highlights the challenge of comparing results from differen t groups using differen t spectra. The 
Li EW estimates for the four cool, main sequence stars in iRandich et al. I (2000) are 8 it 2 mA 
smaller than those found in PRP, leading to an average Li abundance that is 0.07 ±0.01 dex higher 
in the more recent study. Although PRP discounted the memberships of KGP 1405 and KGP 
1365, we include them as possible niemb ers. It also seems probable that PRP followed an incorrect 
identification by INordstrom et al.l (119961 ) of their star 3001 as KGP 1353; it appears that this star 
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i s in f act KGP 1365, with a high probabihty of membership according to iKozhurina-Platais et al 
(jl997l ). 18 stars from PRP overlap with the present study, but only 13 have measurable Li lines 
in both. For these 13 stars, t he equivalent widths in PRP are larger by 4.3 it 11.1 mA, essentially 



midway between the values of iRandich et al. I (|200d ) and PRP, suggesting that our measurements 



are on the same system within the errors. 



Following ISteinhauer &: Delivannis I (|2004l ). we exploit the estimate for the metal abundance 
of the cluster stars to predict and remove the contribution of the Fe I line to the EW measurement 
for Li. The corrected equivalent width for the 6708 A line, listed in Table 2 as EW', along with the 
measured line width, temperature, and s ignal-to-nois e estim ates for each star provide input to a 
grid of abun dance estimates developed bvlSteinhauer I (j2003l ) from MOOG model atmospheres and 
employed by ISteinhauer &: Deliyannis I ()2004l ). For three stars outside the temperature range of the 
model grid, MOOG was used interactively to estimate the Li abundance for the adopted values of 
temperature and surface gravity by comparing synthesized to observed spectra. A measured line 
may produce, at best, an upper lir nit for the Li aburidance if EW' is no larger than three times the 
estimated error, as formulated by iDeliyannis et al. I (jl993l ): stars meeting this criterion are listed 
separately in Table 2. 

A key component in abundance estimation is the scheme by which temperatures are derived 
for stars. For our analysis, temperatures for the unevolved stars have been based exclusively upon 
the redden ing-corrected B — V colors, transformed to the temperature plane using the polynomial 



relation of IDeliyannis. Steinhauer &: Jeffries I (|2002l ): 



Te = 8575 - 5222.7 {B - V)o + 1380.92 {B 



V)o^ + 701.7 {B 



V) ([Fe/H]- [Fe/R]Hyades) 



For our analysis, [Fe/H] for the Hyades was set at +0.15, while the abundance for NGC 3680 
was assumed to be —0.08. For [Fe/H] = 0.0, at a solar Tg of 5770 K, the implied B — V color 
is 0.63, slightly bluer than our normally adopted value of 0.65 for CMD fits to isochrones. If we 
demand that the temperature scale match the redder color, it is equivalent to shifting the scale 
upward by ~65 K. 

Since the relation above is defined by dwarfs in the temperature range bracketing the Sun, it 
isn't applicable to the giants in N GC 3680. For the evolved star s, we converted {B — V)^ to Tg using 
the color-temp erature relatio n of lRa,mirez &: Melendez I ((20051), an update of the dwarf and giant 



calibrations of Alonso et al 



(1996 



d. I1999I1 



For comparison purposes, over the dwarf color range 
of interest, the iRamirez fc: Melendez I (j2005l ) scale predicts temperatures systematically cooler by 
~25 K than the predicted values from the relation adopted above. Final derived values of Tg for 
each star can be found in the second column of Table 2. The derived Li abundances are presented 
in Table 2, with the stars separated according to whether the star had a measurable Li abundance 
or an upper limit as defined by the criterion described earlier. 
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3.3. Li - Comparison with Previous Work 



PR P employed B — F photometry (INordstrom et al.l 119961 ) to estimate temperatm'es for the 



dwarfs (jAlonso et al. 1119961 ) and giants (jAlonso et al. Ill999l ). From the 20 dwarfs and 6 giants 



discussed by PRP, the mean temperature differences, in the sense (PRP - this study), are -83 ± 
63 K and +27 it 57 K, consistent with PRP dwarf colors that are systematically redder by about 
0.03 mag while the giants are bluer by 0.01 mag. Note that the comparison can be made for the 
entire sample because we have B — V colors for every star analyzed by PRP, not just the overlap 
with this work. The offset flip between dwarfs and giants as well as the scatter are products of the 
lower precision B — V colors adopted by PRP. Our colors are based upon a weighted average of a 
number of photoelectric and CCD studies transformed to a common system, resulting in typical 
uncertainties in the photometric colors at the 0.003 to 0.010 mag level. 

To enable a more realistic assessment of the two abundance scales and provide some means of 
merging the two, we have adjusted the Li abundances of PRP to our (B — y)-based temperature 
scale using the prescription for the temperature dependence as supplied by PRP. For seven dwarfs 
with measurable Li in both studies, the mean residual, in the sense (PRP - this study) is -0.02 it 
0.14 dex. If the star with the largest discrepancy is dropped, the average for 6 stars is +0.01 ± 0.10 
dex. The Li abundances for these 7 dwarfs will be based upon the average of the two measures, 
without adjustment for an offset. For the 4 giants with measured Li, the mean residual is -0.11 ± 
0.11 dex. This offset has been applied to the data of PRP and averaged with ours. In the following 
discussions of NGC 3680 as a representative of intermediate-age open clusters, the Li abundances of 
PRP for stars not observed in the current investigation will be used to augment the cluster sample, 
with the offsets based upon the dwarf and giant comparisons applied as needed. 

Based upon multiple abundance reductions for our sample under different assumptions for the 
temperature scale, we find that a change in Tg of ±60 K will alter A(Li) by ±0.05 dex for the dwarf 
stars. For the giants, a comparable temperature change alters A(Li) by ±0.09 dex. 



3.4. Cluster Metallicity - Fe, Ca, Si, Ni 

The determination of the metal abundance of NGC 3680 began with the adoption of a list of 
cleanly measurable iron lines used in common with other cluster investigations by the authors from 
the Indiana University group. To these eight Fe lines we added six measurable lines of neutral Si, 
Ca and Ni. Equivalent widths were measured using the IRAF routine SPLOT, then processed using 
the force-fit abundance routines in MOOG. These steps were replicated for 15 separate exposures 
of the solar spectrum obtained from our daytime sky exposures from the 2005 Hydra run, adopting 
an effective temperature of 5770 K, surface gravity of log g = 4.40 and a microturbulent velocity 
of 1.14 km/sec for the Sun. Signal-to-noise ratios per pixel for each of these solar spectra ranged 
from 90 to 230. All resulting stellar abundances were differenced with respect to the generated 
solar abundances on a line-by-line basis. 
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MOOG provides abundances for each measured equivalent width consistent with an appropri- 
ately selected atmospheric model, so the choices of atmospheric parameters are critically important. 
The choices for temperature scales based on {B — V) colors have been described earlier; gravities 
were determined by examining isochrones fitted to the observed CMD as shown in Fig. 3. For 
composite systems, this will lead to slightly lower gravities than expected for a single star, but the 
change is too small to have a serious impact on the final abundance estimates. Microturbulent ve- 
lo cities for dwarf star s were computed from temperatures and surface gravities using the algorithm 



of lEdvardsson et al. I 



Abundances for the giant stars, for which the lEdvardsson et al. I (|l993l ) algorithm is not ap- 



plicable, are considerably more sensitive to the choice of mi c roturb ulent velocity. As a starting 
point we considered the formula proposed by ICarretta et al. I (j2004l ): vt = 1.5 — 0.13 log g, which 
yields vt values between 1.16 km/sec and 1.24 km/sec for the eight giant stars in NGC 3680. How- 
ever, based upon the standard technique of minimizing the trend between the residuals in [Fe/H] 
and the EW for the lines, PRP found significantly larger values of vt = 1.70 km/sec for their 
coo lest giant and 1 . 90 kin /sec to 1.95 km/sec for the remaining five giants. Attempts to apply 
the ICarretta et al. I (|2004l ) formula to our spectra also failed to achieve appropriate minimization 
of the slope between the residuals and the EW. A range in vt was tested for e ach star until a 



satisfactory slope was achieved. The resulting trend of vt with log g preserved the ICarretta et al. 



(j2004l ) gr avity dependenc e of vt but implied an intercept larger by 0.35 km/sec, midway between 



PRP and ICarretta et al. I (12004 ). 



Results for individual stars are presented in Table 3. For the giant stars, all eight Fe lines and 
the six additional metal lines were measurable in all eight stars. For the dwarf stars, from four to 
six iron lines were measurable in the chosen stars along with the 6717 A line of Ca, the 6721 A Si 
line, and the three Ni lines. The abundances for the dwarf stars exhibited no dependence on Te, 
though the temperature range is modest. An examination of the abundances as a function of signal- 
to-noise, however, suggests that the most consistent dwarf abundances are obtained from the five 
dwarfs with the highest S/N spectra; these five stars are designated in Table 3 by KGP identification 
numbers in bold font. For reference, we note here that we have conformed with what is presumably 
a standard practice when combining separate measures of [A/H], namely constructing a simple 
mean of [A/H] values. It was anticipated that following this practice would facilitate comparison 
with other published work, but it cannot be construed as a correct way to combine logarithmic 
values. The preferred approach (taking the logarithm of averaged numerical abundances) yields 
abundances that are typically larger by amounts that reflect the scatter among individual values. 



3.5. Metallicity - Giants 

The mean [Fe/H] for 8 giants in Table 3 is -0.17 it 0.07, with comparable values for Si of +0.05 
lb 0.12, based on one line, and -0.17 it 0.08 for the three lines of Ni. A weighted average for [Fe/H] 
gives -0.18 lb 0.07. While the balance of the measures clearly favors a subsolar abundance for NGC 
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3680, it leads to the obvious question of the zero-point of the metalhcity scale since the solar spectra 
and their analyses are distinctly removed from the much cooler and more luminous giants. To test 
our [Fe/H] estimate, we turn first to the six giants analyzed by PRP, five of which have abundances 
in our survey. PRP obtained [Fe/H] = -0.26 it 0.11 from 6 giants. The first step in the comparison 
is to correct for the differences in temperature caused by PRP's use of less precise photometry, a 
smaller reddening, and a different B — V calibration. This adjustment produces only a modest 
decrease in the PRP mean value to [Fe/H] = -0.28, while reducing the scatter to =b 0.08 dex. The 
truly significant adjustment comes from the difference in the adopted vt- For the 5 stars common 
to the two samples, if we had adopted the microturbulence values used by PRP, the mean [Fe/H] 
would shift by -0.13 ± 0.05 dex. This is exactly the difference and the scatter one gets by taking 
the average residual for the 5 giants after placing the abundances on the same temperature scale. 
We have shifted PRP's data for the 6 giants by +0.13 dex and averaged our values with theirs for 
the 5 stars in common to derive a mean [Fe/H] = -0.15 it 0.08 from 9 giants. The slightly higher 
value of [Fe/H] comes from the inclusion of KGP 1379 (E34), which has the highest abundance 
of all the giants. We opted not to analyze this star since its composite color would likely distort 
any temperature estimate. As one might expect, if this star is dropped, the mean [Fe/H] of the 
remaining 8 giants becomes -0.17 ± 0.06. 

Two high-resolution s p ectro scopic investigations of giants in NGC 3680 have been published 



recently. 
[Fe/H] = 



Smiljanic et al. 



(120081 ^ discussed the giant KGP 1175 (EG13), for which they derive 
+0.04 ±0.11 (std dev). In spite of their note to the contrary, this star is one of the five 
studied by PRP that overlaps with the current investigation. For their analysis, the stellar param- 
eters Tg = 4660 K, log g = 2.6, and vt = 1.3 km/sec were derived purely from the requirements of 
internal consistency among the line abundances with EW, ionization state, and excitation poten- 
tial. Our temperature for this star is 4606 K, while PRP get 4668 K. Adjusting to our scale would 
lower the [Fe/H] to approximately -0.01. Our log g is 2.56, so no correction is required, leaving a 
difference of 0.12 dex between their estimate and ours for this star. If our derived vt of 1.56 km/sec 
is reduced to 1.30 km/sec, our [Fe/H] is shifted by 0.11 dex, accounting for the remainder of the 
difference. 



The second paper (ISantos et al. II2009I ) details spectroscopy of three single-star giants, 1175, 
1374, and 1461, analyzed using two different line lists, with the stellar parameters set spectroscopi- 
cally. The two line lists generate different combinations of the stellar parameters for the same star, 
but the systematic offset between the abundances for the three giants is only 0.02 dex, even though 
the temperature scales differ by 127 K. The temperature shift appears to compensate in part for 
the systematic change in vt, which drops by an average of 0.05 km/sec for the lower Tg scale. As 
an ex ample of the s i ze of the effects we are discussing, for KGP 1175 on the cooler temperature 
scale. ISantos et al. I (|2009l l find (4657 K, 2.68 dex, 1.3 8 km/sec) and [Fe/H] = -0.01. The expected 



value with adoption of the vt lowered to 1.30 km/sec (ISmiljanic et al. 1120081 ) is +0.03. 



If we adopt the parameters tied to either temperature scale of ISantos et al. I (|2009l ) for ana- 



lyzing our spectra for the three giants common to the two studies, [Fe/H] estimates for KGP 1175, 
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1374, and 1461 are increased by +0.12 it 0.02 dex for the cooler scale and +0.15 ± 0.02 dex for the 
hotter scale. If the offsets are applicable to the entire giant sample, our met allicity from th e giant s 



becomes [Fe/H] = -0.05 and -0.02, respectively, the same values found by ISantos et al. I (|2009l ) 



Note that the very good agreement among the abundance estimates when the same stellar param- 
eters are applied weakens the argument that the discrepancy among the studies is a byproduct 
of underestiniated E W from spectra of lower resolution and S/N compared to those available to 



Santos et al. I ([20091). 



The good news, therefore, is that the high-dispersion spectra for giants discussed to date, when 
processed using the same fundamental parameters, will generate Fe abundances that are similar 
at a level typically ±0.02 dex. The bad news is that for a fixed temperature scale, the absolute 
abundance scale is dominated by the adopted microturbulence velocity, vt, which leads to a range 
from [Fe/H] = 0.05 to —0.30 for vt between 1.3 to 1.9 km/sec. 

A potential handle on the location of the zero-point of the metallicity scale comes from the 
observation by PRP of two stars in the Hyades, for which they der i ve [Fe /H] = +0.06. Since the 
temper ature scale for the Hyades gi ants as defined bv lAlonso et al. I (|l999l ) is virtually identical to 
that of iRamirez &: Melendez I (|2005l ). the modification applied to the NGC 3680 abundances due to 
the altered photometry, reddening, and metallicity doesn't apply. Thus, any remaining offset should 
be a byproduct of the adopted vt- Assuming the same formulation applied to derive vt for NGC 
3680 results in the same systematic offset when applied to the Hyades, if one adopts a metallicity 
for the Hyades near the higher end of the traditional values, [Fe/H] = +0.15, this shifts the PRP 
corrected abundance for NGC 3680 from -0.29 to -0.20, still out of range of solar abundance but 
a respectable match to the value derived in this study. We will return to this question again after 
the discussion of the results for the dwarfs. 



3.6. The Dwarfs 



For the 13 stars at the turnoff with S/N above 70, the average [Fe/H] = -0.04 ± 0.11, while Si, 
Ca, and Ni have [A/H] = -0.06 ± 0.10, -0.17 ± 0.12, and -0.12 ± 0.11, respectively. As noted earlier, 
the scatter among the abundances is predictably correlated with S/N. If the analysis is limited to 
the 5 stars in Table 3 with S/N above 110 (boldface ID), the Fe, Si, Ca, and Ni abundances become 
-0.07 ± 0.02, -0.13 ± 0.05, -0.22 ± 0.07, and -0.18 ± 0.11, respectively. Within the uncertainties, 
the abundances are consistent with a scaled-solar distribution. 



For comparison. iPace et al. I (|2008l ) and I Santos et al. I (j2009l ) have derived abundances for two 
G-dwarfs in the cluster from the same spectra, finding [Fe/H] = 0.00 and -0.03 for KGP 1752 and 
-0.07 and -0.04 for KGP 1050, respectively, leading to an average of [Fe/H] = -0.04, consistent with 
either value for our turnoff dwarfs. However, for consistency, these estimates deserve closer exam- 
ination. Unlike the giants, the errors in the dwar f abundance s are d ominated by the temperature 
scale and are only very weakly dependent upon vt- iPace et al. I (|2008l ) use three reddening-corrected 
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color indices coupling to 2 MASS JHK photometry, then transferred to Tg using the relations of 
Ramirez &: Melendez I (j2005l ) and averaged to define a photometric temperature. This value sup- 
plies the starting point for a grid of stellar atmospheres that are interpolated among the variables 
until the line abundances show minimal trends with EW, excitation potential, or ionization state. 



For star KGP 1050, iPace et al. I (120081 ) determined the photometric temperature to be 6053 
K, similar to our B — V based estimate of 6098 K, but hotter than the value adopted in PRP by 
174 K. However, the spectroscopic temperature used to estimate the abundance was 6210 K, 112 
K hotter than the scale used in our abundance derivation for the dwarfs. Use of our photometric 
scale would lower the abundance of KGP 1050 to [Fe/H] = - 0.16; using the photo metric scale of 
Pace et al. I (|2008l ) would shift the value to -0.20. By contrast, ISantos et alTI ^QM) adopt 6134 K, 
36 K hotter than our photometric value; use of our scale would reduce their abundance to [Fe/H] 
= -0.07. 

Star KGP 17 52 follows a similar pattern in that the spectroscopic temperature (6010 K) 
(jPace et al. II2008I ) is 84 K higher than the VJHK average (5926 K), but in this case the photo- 
metric value i s only 10 K hotter than the PRP temperature and 6 K cooler than the estimate in 
Santos et al. I (|2009l ). We have no direct means to tie the star reliably into our temperature scale 
because its composite B — V color has a dispersion above 0.05 mag, one of the largest among the 
brighter stars in the photometric catalog. If we assume that the VHK va lue is compa r able t o our 
system, the abundan ce for this star would decline to [Fe/H] = -0.06 for iPace et al. I ()2008l ) and 
remain unchanged for lSantos et al. I (|2009l ). 



However, there is a more important reason why the results for this star remain questionable 
and why it wasn't included in our investigation. The star has a high proper-motion membership 
probability of 58% a nd its radial velo c ity is consistent with being a single-star member of the 
cluster. Despite this, iNordstrom et al.l (119971 ) flag the star as a questionable member due to its 
suspicious position in the C MP. The photometric and spectroscopic temperatures for KGP 1050 
and 1752 (jPace et al. 1120081 ) imply a temperature dif ference between thes e two stars of 127 K and 
100 K, respectively, while the spectroscopic scale of Santos et al. (j2009l ) leads to a difference of 
202 K. A typical shift of 110 (200) K requires KGP 1752 to be 0.03 (0.05) mag redder than KGP 
1050. On the main sequence, being cooler by this amount should place KGP 1752 almost 0.2 (0.3) 
mag fainter than KGP 1050. Instead, KGP 1752 is 0.3 mag brighter than KGP 1050, placing it 
0.5 (0.6) mag above its expected position in the CMD for a star of its color. If the star is not a 
binary (for which there is no spectroscopic evidence), it is either a non- member or its structure does 
not resemble that of a typical solar-type star. For either of the above reasons, any spectroscopic 
abundance tied to this star remains questionable. 

What is the correct metallicity for NGC 3680? From the data obtained solely from this study, 
the 8 giants have [Fe/H] = -0.17 it 0.07 while the best data for the 5 high S/N dwarfs imply [Fe/H] 
= -0.07 lb 0.02, with a simple average of -0.12. If one derives a weighted average tied to the inverse 
of the standard errors of the mean, [Fe/H] = -0.09. When analyzed adopting the same parameters 
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selected in this investigation, the 6 red giants of PRP and the one giant of ISmilianic et al. I (|2008l ) 
have [Fe/H] = -0.15 ± 0.08 and -0.13 ± 0.10, respectively. As emphasized previously, the problem 
with the giants is that by an appropriate choice of vt , the mean metallicity could easily be raised to 
solar or higher. The one definitive dwarf member an alyzed outside this investigation, when placed 
on our temperature scale, has [Fe/H] between -0.07 (jSantos et al. II2009I ) and -0.16. 



From past work, the excl usive me tallicity estimate for the dwarfs for over 15 years was the 
uv byil0-hased value of -1-0.1 of iNissen reproduced u sing the same photometric zero-points 



bv lAnthonv-Twarog et al. I (|l989l l and Bruntt et al. I (|l999l l. More recently, the uvbyCaRP CCD 



study of ATT identified a possible error in the mi photometry of the earlier work and obtained a sec- 
ond metallicity estimate using the independently calibrated hk index. From 34 stars at the turnoff 
using both photometric indices, mi and hk, the cluster metallicity is found to be [Fe/H] = -0.14 
lb 0.03 on a scale where M 67 and the Hyad es are found to have [ Fe/H] = -0.06 and -1-0.12, respec- 
tively (INissen et al. 1119871 ). For the giants. iTwarog et al. I (Il997l'l used the combi ned observations 
on the DDO system with the moderate dispersion spectra of lFriel Sz Janes (1993), transformed to 
a common metallicity scale, to obtain [Fe/H] = -0.10 it 0.06 from 12 measurements. On the same 
metallicity scale tied to DDO photometry an d moderate-dispersion spectra, 36 measurements for 
giants in M67 produce [Fe/H] = 0.00 ± 0.09. iTwarog et al. I (jl997l ) also provide a possible avenue 



for constraining the zero-point sca. t ter am ong the various high dispersion abundances for the cluster 
dwarfs and giants. ISantos et al. I (|2009l ) investigate the relative abundances between dwarfs and 
giants in 13 open clusters, finding a small effect with met a llicity that is readily corrected. Of the 13 
clusters, 12 are included in the catalog of iTwarog et al. I (|l997l ). For these 12 clusters, the average 
difference in [Fe/H], in the sense (santos - taat), is 0.033 it 0.050 dex. The expected scatter from the 
quoted errors for the individual abundance estimates is ±0.046 dex. If the residuals are weighted 
by the inverse erro rs, the mean of f set ris es to 0.042 dex. With [Fe/H] = -0.04 from 3 giants in NGC 



3680, the result of ISantos et al. I (j2009l ) becomes [Fe/H] = -0.08 on the system of ITwarog et al 
(|l997l ) where, again, the NGC 3680 giants were found to have [Fe/H] = -0.10 ± 0.02 (sem) from 
12 spectroscopic and DDO measurements. 

An average of the various abundance estimates without regard to zero-point and weighted by 
the inverse of the errors is dominated by the dwarfs: [Fe/H] = -0.08 it 0.02 (sem). Adoption of 
values between [Fe/H] = -0.12 and -0.04 have little impact on our conclusions regarding the Li 
abundances. A change in [Fe/H] of ±0.04 will have a small effect on the ages (±0.1 Gyr), distance 
moduli (±0.05 mag), and Li-dip masses (±0.03 Mq). 



4. The Li-Dip: Cluster Data 
4.1. NGC 3680 



With the establishment of a plausible set of cluster properties for NGC 3680, we can now 
attempt to place the cluster's Li-dip in the context of the physical parameters that define the dip 
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and its evolution ov er time. The key pape rs that have dealt with th is issue previ o usly i n a compre- 



hensive manner are iBalachandran I (|l995l ) using open clusters and IChen et al. I (120011) eva, l uating 



field dwarfs. In a discussion of the pattern of A(Li) in several open clusters, IBalachandran I (|l995l ) 
concludes that the masses of the stars in the dip depend on metallicity, but the ZAMS surface tem- 
peratures of the stars within the dip do not. The temperature range of the dip appears to depend 
less on metallicity than on age, based on a comparison of Praesepe and the Hyades, two clusters 



of sim ilar age but supposedly different metallicity, a point we will return to below. IChen et al. 



(|200ll ) similarly find that the location of the dip is metallicity-dependent and derive a relation for 
the mass of the central dip stars as a function of [Fe/H], with solar metallicity dip stars having 
masses near 1.4 Mq, declining to 1.1 Mq at [Fe/H] = -1.0. 

The canonical illustration of the main sequence Li-dip is a plot of Li abundances versus main- 
sequence surface temperature or color. For clusters as old as NGC 3680, beyond illustrating the 
difference in Li between the main sequence and the evolved giants, this representation provides 
little insight since evolution has shifted the Li-dip into the vertical strip at the CMD turnoff 
(jPaniel et al. I Il994l : IBalachandran I Il995l ) . The more functional approach is to use V as the dis- 
criminator of mass, keeping in mind the potential confusion caused by binary systems. The relevant 
diagram using our merged and homogeneous sample for NGC 3680 is shown in Fig. 4. Filled cir- 
cles and triangles represent stars with detectable Li and upper limits, respectively, for stars in our 
sample, averaged with PRP where available. Open circles and triangles show stars with detectable 
Li and upper limits, respectively, from the sample of PRP alone. With the sample boosted from 18 
to 40 stars, the features defining the Li distribution among the turnoff stars in NGC 3680 are en- 
hanced and expanded. Keeping in mind that binarity can shift the position of a star in the diagram 
by, at most, 0.75 mag to the right and the predominance of upper limits for Li estimates within 
the dip region, the location of the center of the dip is reasonably defined as ^ = 13.45 it 0.15. 
This approximate location is obtained by giving significant weight to the distribution of measured 
Li abundances that comprise the wings of the dip rather than the upper limits that dominate the 
core of the dip. The midway point between the wings is between V 13.5 and 13.6. Taking into 
consideration the fact that the drop in Li is normally somewhat steeper on the high-mass side, this 
was adjusted to a central dip at 13.45, though clearly values as high as 13.6 are plausible, as defined 
by the error estimate. A more exact definition of the location of the dip center would require an 
increase in cluster members (not possible) or the superposition of data from clusters of comparable 
age and composition, the focus of Section 5. 

From the minimum in the Li-dip, the Li abundances for higher mass stars rise steeply to an 
upper bound of A(Li) = 3.35 hy V = 12.85, maintaining this high value to at least V = 12.4. There 
is one star (KGP 1035) at the Li upper bound at V = 11.9; it is classed bv lNordstrom et al.l (|l997l ) 



as an SB2 and its position in the CMD places it in an evolutionary phase that is so rapid that it is 
highly likely this star has been shifted to smaller V from a true position between V = 12.4 and 12.8 
by its composite nature. By contrast, the one star brighter than the red hook in the turnoff which 
shows no evidence of binarity (KGP 1522) has A(Li) = 2.80, almost a factor of 4 below the upper 
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bound. This is consistent with a significant Li reduction as stars evolve through the hydrogen- 
exhaustion-phase toward the red giant branch, where t he highest Li abundan ces are reduced by 



over two orders of magnitude relative to the upper bound (jPasquini et al. II2004I ). When contrasting 
observation and theory for clusters of this age, it's important to remember that the exact shape 
of the turnoff hook depends very sensitively on the treatment of convective core overshoot in the 
models. 

On the lower mass side of the dip, the Li abundance rises almost as rapidly, reaching a maximum 
at A(Li) = 2.75 it 0.15 near V = 14.2 it 0.1, with seemingly little scatter. Toward fainter V, the 
Li abundance declines by an order of magnitude, reaching a value near A(Li) = 1.8 hj V ^ 15.2. 

For purposes of constraining the dip and its evolution, the V dependence of the 1.8 Gyr 
isochrone in Fig. 3 has been converted to a mass scale, as illustrated along the top of the plot in 
Fig. 4. Based on the V-mass relation, the ZAMS mass of the stars that occupy the center of the 
Li-dip is 1.35 it 0.03 Mq. The current edges of the dip range between 1.51 it 0.03 Mq at the high 
end to 1.19 it 0.02 Mq at the low end. If we were to transform this mass range to a main sequence 
at [Fe/H] = -0.08 at the age of the Hyades (0.7 Gyr), the colors coupled to these mass points would 
he B-V = 0.40, 0.32, and 0.50, respectively. 



4.2. NGC 752 



The first and most straightforward comparison for NGC 3680 is with the very similar open 
cluster, NGC 752. Building upon the extensive discussion in ATT, the key issue in this, as in 
any, comparison is how similar. There is little disagreement over the relative reddening of the two 
clusters. Typical estimates for NGC 752 range between E{B - V) = 0.02 and 0.04, while NGC 
3680 currently is bracketed by 0.05 to 0.07. We wih assume E{B - V) for NGC 752 = 0.035, 
consistent with our earlier work. The more contentious issue is the metallicity. An early evaluation 
of the state of the metallicity measures for the clus ter includ i ng st andardization and reanalysis 
of all available photometric data was undertake n bv iTwarog I (jl983l ) and updated and expanded 



to include spectroscopic data by iDaniel et al. I (119941 ). leading to a composite [Fe/H] of -0.15 



with significant weight given to the spectroscopic results. Photometric zero-point issues aside, this 
conclusion would imply an abundance comparable to NGC 3680. Due to improved calibrations, 
better zero-point photometric checks, and expanded spectroscopic data, the earlier results have 
gradually been superceded, with a modest change in the answer. Using the co mbined homogeneous 



results from DDO photometry and moderate-dispersion spectra of the giants (jTwarog et al. 1119971 ) 



in the manner outlined in Sec. 3.6 for NGC 3680, the metallicity of NGC 752 is found to be [Fe/H] 
= -0.09 lb 0.07 from 16 measures. It should be emphasized in this particular case that, unlike 
NGC 3680, the DDO result (-0.16) differs from that of the spectroscopy (-0.05) by more than 0.1 
de x. A sigii i ficant improvement occurs in the results from uvbyllf3 ph otometry, where the m i data 



and 



of iTwarog I (119831) exhibited a systematic offset compared to that of ICrawford &: Barnes I (jl970l ) 



Nissen I (jl988l ). thereby generating an abundance below [Fe/H] = -0.2. With the resolution 
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of the zero-point question (jjoner &: Taylor Ill995l ). the mean [Fe/H] of NGC 752 on the uvbyHlS 
system is well established at -0.08 it 0.11 from 31 stars, on a scale where NGC 3680 has [Fe/H] = 
-0.17 from mi data alone and M67 and the Hyades are at -0.06 and +0.12, respectively. Recently, 
Anthony- Twarog &: Twarog I (|2006l ) observed NGC 752 for the first time on the extended uvbyCa 
system and, from 10 single-star members at the cluster turnoff, find [Fe/H] = -0.06 it 0.09 on a 
scale where the Hyades has [Fe/H] = +0.12. 

Despite its proximity to the Sun, there have been surprisingl y few spectroscopic studies of stars 
within NGC 752. At high dispersion. [Hobbs &: Thorburn I (|l992l ) find [Fe/H] = -0.09 ± 0.05 from 8 
turnoff stars on a scale whe re similar analysis in M67 generates [Fe/H] = -0.04 ± 0.12 for 5 dwarfs. 
The high-dispersion data of lSestito et al. I (|2004l ) for three stars imply [Fe/H] = +0.01 ± 0.04 from 
7 dwarfs, on a scale where th e Hyades has [Fe/H] = +0.14 ± 0.06, based on 3 stars. A similar 



analysis from the same group (IRandich et al. II2006I ) leads to [Fe/H] = +0.03 ± 0.03 froni 10 stars 



in M67 on a scale where one Hyades star gives [Fe/H] = +0.13. While ISestito et al. I (120041 ) dismiss 
all previous metallicity estimates as either irrelevant or unreliable, a careful analysis confirms that 
they are remarkably consistent, with most comparisons of the absolute scale distorted by differences 
ai nong the metal li city z ero-p oints. A classic ex ample is presented by the moderate-dispersion work 
of iFriel &: Janes I (|l993l ) and lFriel et al. I (|2002l ) where cluster abundances are derived with a zero- 
point for the metallicity scale defined so that M67 has [Fe/H] = -0.09 in the original analysis and 
is fixed at [Fe/H] = -0.15 in the latter paper. Any comparison of the results from these papers to 
a spectroscopic scale where M67 is observed to have [Fe/H] = +0.03 will lead automatically to a 
spurious discrepancy if no zero-point correction is applied. 

From the discussion above, the qualitative pattern that emerges is that NGC 752 is more 
metal-rich than NGC 3680, but more metal-poor than M67. A weighted average of the above 
techniques leads to [Fe/H] = -0.05 ± 0.04 on a scale where M67 has [Fe/H] = 0.00 and the Hyades 
is at +0.12. 

With the reddening and metallicity differential set, we can now make use of the results of ATT 
to place the clusters on the same apparent magnitude scale. By direct superposition of the cluster 
CMD's, ATT found that the apparent modulus of NGC 752 had to be 1.90 mag smaller than that 
of NGC 3680, assuming they had the same metallicity. With a higher [Fe/H], the main sequence 
of NGC 752 should be brighter than that of NGC 3680 by approximately 0.03 mag, reducing the 
offset to 1.87 mag. With the derived (m - M) for NGC 3680 set at 10.30, this imphes {m - M) = 
8.43, a boost from the value of 8.30 found in ATT. Using an appropriately interpolated set of 
isochrones for the slightly higher metallicity adopted for NGC 752, the age of the cluster becomes 
1.45 ± 0.1 Gyr, as shown in Fig. 5, a slight reduction from the value of 1.55 Gyr derived in ATT. 
Note that the typical position of the giants is blueward of the first-ascent relation by 0.07 mag. 

The combined Li data for NGC 3680 (black symbols) and NGC 752 (red symbols) are shown 
in Fig. 6, with the data for NGC 752 shif t ed to t he distance of NGC 3680. Th e dat a for NGC 752 



are compiled from iHobbs &: Pilachowski I (jl986l ) , iPilachowski &: Hobbs I (jl988l ) and ISestito et al. 
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(j2004l ). Because each of these NGC 752 studies hsted equivalent widths for the Li hne, we were 
able to analyze the data in a manner consistent with our treatment of our own measures of Li in 
NGC 3680, including use of the temperature-color calibration described above. A metal abundance 
[Fe/H] = -0.05 was assumed f or this purpose, with a reddening value of E{B — V) = 0.035 applied 
to {B — V) colors compi led bvlDaniel et al. I (| 19941 ) . For two of the t hree stars for which equ i valent 
widt h measurements by 



Sestito et al. 



(j2004l ) replicated measures by iHobbs Sz Pilachowski I (jl986l ) 
and iPilachowski &: Hobbs I (Il988l ). the resulting Li abundance estimates were averaged. For one 
of the stars, P859/H207, the published equ i valen t widths are widely discrepant. The measured 
equivalent width published by ISestito et al. I (|200J) was used to compute the Li abundance for this 
star. 

From Fig. 6, the morphological agreement between the two clusters is impressive. There is a 
clear gradient in the upper bound relation as one moves from the brighter stars to the fainter main 
sequence. The Li abundance pattern near the solar system Li value is, within the errors, the same 
for both clusters among the higher mass stars above the gap. The sample of stars in NGC 752 with 
measured Li in this CMD region is frustratingly small; we would predict that a larger sample of 
observations of stars at the turnoff for V brighter than 11.4 in NGC 752 would produce a pattern 
identical to that in NGC 3680 unless parameters other than age and metallicity are at work. 

Moving down the turnoff, the Li-dip is defined at the high-mass edge by a steep drop at V 
= 13.2, in contrast with the low-mass edge which shows a gradual rise from V = 13.6 to A(Li) = 
2.8 at V = 14.25. Toward lower masses, the small sample in NGC 3680 meshes perfectly with the 
decline observed in NGC 752. One potential point of concern could be the fact that the Li-dips 
superpose so smoothly even though NGC 3680 is 0.3 Gyr older than NGC 752 on the current scale. 
However, the mass of the stars within the Li-dip is low enough that between 1.45 Gyr and 1.75 
Gyr the position of the dip shifts to brighter My by less than 0.05 mag. Taking into account the 
slightly higher metallicity of NGC 752 relative to NGC 3680, the mass of the dip stars in NGC 752 
should be slightly higher, so the stars should be slightly more evolved at the same age compared 
to a cluster with lower [Fe/H]. Thus, we expect that the metallicity effect decreases the impact of 
the age difference to a level which is below detectability for the current sample. 



4.3. IC 4651 

As noted in the Introduction, any discussion of NGC 3680 or NGC 752 invariably leads back 
to the third member of the trio, IC 4651. While a common approach is to superpose the Li data 
of IC 4651 upon that of NGC 3680 and NGC 752, the higher metallicity of this cluster provides 
an equally important handle on the evolution of the Li-dip and the reliability of predictions based 
upon the adopted theoretical isochrones. Based on the masses and projected evolution of stars in 
the Hyades Li-dip and turnoff, are the observations of the Li-dip stars in the ~1 Gyr older but 
chemically similar IC 4651 consistent? 
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The preliminary step is to place IC 4651 on the same parametric scale as NGC 3680 and 
NGC 752. Unlike NGC 3680, there has been remarkable consistency in the estimation of the 
cluster parameters from both photometry and spectros copy. On the photometric side, irideperi - 
den t investigations of the clus t er ori t he uvbyil0 system. lAnthonv-Twarog &: Twarog I (jl987l . [2000) 



and iNissen I (|l988l ) ; iMeibom I (j2000l ) ; iMeibom et al. I (j2002l ) , have exhibited none of the zero-point 
issues identified for NGC 3680 and have always indicated that the cluster is similar in metallic- 
ity to the Hyades. Equally consistent ar e the reddeni n g est imates tied to the H/3 photometry of 
Anthony- Twarog &: Twarog I (|l987l . 12000) and iNissen I (119881 ) . When analyzed using the same in- 
trinsic color relations, the two data sets produce E{b — y) = 0.071 and 0.076, respectively, equiv alent 
to E(B — V) = . 10 an d 0.11. With these reddening values. I Anthonv- Twarog fc Twarog I (|2000l ) and 
Meibom et al. I (120021') find [Fe/ H] = -1-0.12 and and -1-0.13, respectively, the latter after adjusting 



the nil values of INissen I (jl988l ) downward by 0.005 mag to coincide with the their CCD survey. 
Likewis e, DDO data for 11 g iants leads to E{B — V) = 0.11, and a corresponding [Fe/H] = -1-0.10 



± 0.05 (|Twarog et al. Ill997l ^. 



On the spectroscopic side, ICarretta et al. I (120041 ) analyzed 4 giants to derive [F e/H] = -1-0.11 ± 
0.01, a ssuming E{B—V) = 0.083; adoption of a higher reddening raises this estimate. iPasquini et al. 



(j2004l ) discuss high-dispersion spectroscopic analysis of over 20 stars ranging from giants to G dwarfs 
and obtain [Fe/H] = +0 10 ± 0.03 with an ad opted E{B - V) between 0.12 and 0.13. Finally, from 



(120091) derive [Fe/H] = -^0.15 ± 0.01 and +0.15 + 0.02, 



3 giants and 5 dwarfs, 1 Santos et al. 
respectively. 

For the apparent distance modulus, the consistency among recent determinations tied to main 
sequenc e fitting has be e n equ ally impressive, with one notable exception which will be discussed 



below. iTwarog et al. 1 (|l997l ) find (m - M) = 10.25 for [Fe/H] = +0.1 and E{B - V) = 0.11 
and ATT derive (m — M) = 10.20 with reddening reduced to 0.10. Given the similarity of the 
cluster metallicity to the Hyades, one can avoid the usual issues with isochrqnes b y doing a direct 
comparison to the Hyades. Using the Vby data. lAnthonv-Twarog &: Twarog I ({2000 ) find (m - M) 
= 10.30 for E(b — y) = 0.(3 71 {E{B — V) = 0.10) and from the Vvy diagram to much fainter 
magnitudes. IMeibom I (|20od) finds {m - M) = 10.36 for E{b - y) = 0.076 {E{B -V) = 0.11). 



While the cluster parameters are exceptionally consistent internally, one minor discordant 



note has been raised by rec ent high-dispersion spe ctroscopic work on the cluster. iPasquini et al. 



(l2004l ). lBiazzo et al. 1 (120071 ) and lPace et al. I (120081 ) have used temperatures based upon consistency 
constraints on abundances from the multiple lines of differing equivalent width, ionization state, 
and excitation potential to imply that the reddening value of IC 4651 is closer to E[B — V) = 
0.12 to 0.13. This co nclusion is contradicted by a similar ap proach applied to the red giants by 
Carretta et al. 1 (1200J) and severely weakened by the analysis of lSantos et al. 1 (l2009l ) , who generate 
the same abundances but two distinctly different temperature scales for the same stars using two 
different line list s. Ultimately, the differences between the higher reddening value and E{B — V) = 
0.10 discussed in lAnthonv-Twarog fc Twarog I (I2OO 0'. '2OO4I) : ll^ibom et al. \hooi \ and E{B-V) 



0.11 discussed in Meibom 



( 2OO0I ) and Twarog et al. 1 ((1993) are minor and well within the error bars 
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of both approaches. If we adopt J^{B — V) = 0.12 for IC 4651, E(6 — y) becomes 0.085; shifting from 
the older average of E(6 — y) = 0.073 leads to a correction to mi of an additional 0.004 mag, raising 
[Fe/H] for IC 4651 on the Stromgren system to [Fe/H] = +0.16, again, not a statistically significant 
change from earlier results. With a reddening of E(i? — V) = 0.12, the apparent modulus rises to 
{m — M) = 10.40 lb 0.05. Under the assumption that IC 4651 has a metallicity indistinguishable 
from that of the Hyades and using an appropriately interpolated set of isochrones ([Fe/H] = 
0.15), E(B -V) = 0.12 and (m - M) = 10.40, as illustrated in Fig. 7, the age of the cluster 
becomes 1.50 j = 0.1 G yr, a reduction from the value of 1.75 Gyr derived in ATT and 1.7 Gyr in 
(|2002l ) due to the higher adopted reddening, a slightly higher absolute [Fe/H], and 



Meibom et al. 



exact interpolation of the isochrones. Triangles are stars known to be spectroscopic binaries; known 
non-members ha ve been removed from the samp le which is a combination of the BV photographic 
and CCD data in I Anthonv-Twarog et al. I (|l988l ). It should be noted that the color and magnitude 
offsets applied to the isochrones for NGC 3680 have also been applied to the se ts at higher [Fe/H] . 
A consistency check has been made by matching the photometry of the Hyades (|Joner et al. 1120061 ) 
to an isochrone of appropriate age and {m — M) = 3.33, with the offsets included; the agreement 
is excellent. 

The comparison between the data and the isochrones is encouraging, especially at the red hook 
in the turnoff which is dominated by the assumed degree of convective overshoot mixing. Note, 
however, that the giants all lie blueward of the first-ascent giant branch, in sharp contrast with NGC 
3680 and as predicted if the giants are dominated by He-core-burning stars. This difference with 
NGC 3680 could be removed by altering the reddening of IC 4651 to a value closer to E(i? — V) = 
0.00, thereby requiring a reduction in the apparent modulus by 0.6 mag and significantly boosting 
the cluster age. As discussed previously, however, reductions in the reddening go opposite to the 
trend from spectroscopy in recent years. 

Before the IC 46 51 data are lin k ed to younger clusters of the same metallicity, the particularly 
anomalous results of iBiazzo et al. I (j2007l ) need to be addressed. The focus of their study is the 
use of line-depth-ratios from high-dispersion spectra of giants to derive temperatures which, when 
combined with absolute magnitude and metallicity, car i generate a ru a ss an d, ultimately, through 



isochrone comparisons, an age. The problem is that iBiazzo et al. I (|2007l ) quote an age of 1.4 
Gyr and an apparent modulus of (m — M) of 9.83 for IC 4651 from a comparison to the Padova 
isochrones. We have reproduced the comparison of Fig. 7 using properly interpolated and zeroed 
Padova isochrones and find excellent agreement with the age a ,nd modulus o b taine d from the Y"^ 
isochrones, as expected from the comparison with ATT. The IBiazzo et al. I (|2007l ) age estimate 
should be adjusted downward to correct for their use of [Fe/H] = 0.10, if their results are to be 
compared with our adopted [Fe/H] = 0.15. Clearly, however, an estimate of (m — M) = 9.8 is 
implausible given the discussion above and the fit in Fig. 7; th e minimum modulus using the lowe r 
limit to the redd ening, 'E(B — V) = 0.08 5, is 10.15, as derived in I Anthony- Twarog &: Twaroe I (120001) 



and a dopted by ICarretta et al. I (j200J) , and can only increase for higher reddening. 



(|2007l ) make no note of this apparent anomaly, citing only two dated references (jNissen 



Biazzo et al. 



1988 
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Kieldsen Frandsen 1 ll99lh that are consis tent with thei r resul t. The primary source of the problem 



comes from their use of the photometry of iPiatti et al. I (119981) to compare the cluster CMD to the 



theoretical isochrones. As discussed in lAnthonv-Twarog &: Twarog I (|200Cl ). the V magnitudes in 
this survey are systematically too bright by 0.20 mag. Even worse, the V zero-points between the 
two CCD fields studied in the cluster differ by 0.09 mag. Because no additional V — I colors are 
available for comparison, any attempt to use their V — I CMD to derive an age by matching to 
isochrones must be highly questionable, making the approximate agreement between 1.4 Gyr and 
Fig. 7 purely fortuitous. 



4.4. Hyades and Praesepe 



As a younger Li standard for comparison purposes, we have combined data for the Hyades 
and Praesepe. There is no question that the Hyades is metal-rich, with [Fe/H] values noted in this 
paper between -1-0.12 and -1-0.17; we have selected -1-0.15 for the temperature scale. What matters, 
however, is that IC 4651 and the Hyades be the same within ±0.05 dex. The more controversial 
cluster is Praesepe. Following the same sources discussed for the o ther clusters, DDO p hotometry 



of 4 giants pro duces [F e /H] = +0.U ± 0.07 and E{B - V) = 0.0 (jTwarog et al. II1997I ). uvbyRp 
photometry of iNissen I (jl988l ) implies E{B - V) less than 0.01 and [Fe/H] = +0.W ± 0.15 from 
42 dwarfs on a scale where the Hyades has [Fe/H] = -1-0 .12, consistent with the earlier results of 



index of 



Crawford &: Barnes 



Crawford &: Barnes I (119691). iJoner &: Tavlor I (|l995l . 120071 ) derive an offset of ~-|-0.01 mag to the mi 



(jl969l ). which would imply a higher [Fe/H]. Washington photometry 
of 4 giants in each of the Hyades and Praesepe clusters leads to +0.07 and -1-0.04, respectively 



(|Geisler et al. Ill99lh 



Sp ectroscopically, the high-dispersion staridard f or some time has been t ied to three related 
papers (jBoesgaard fc Budge I Il988l : iBoesgaard I Il989l : iFriel Boesgaard I Il992l ) . The first investi- 
gation found [Fe/H] = -1-0.13 =b 0.07 (sem) from 5 single dwarfs, on a scale where F dwarfs in the 
Hyades had [Fe/H] = -|-0.17 =b 0.06 (sem). The latter two studies generate a combined [Fe/H] = 
-1-0.06 =b 0.06 from 7 dwarfs on a scale where the Hyades has [Fe/H] = -1-0.13; a weighted averag e 



combining the se two saraples r aises the Praesepe value to -1-0.09 ± 0.03 (sem) (lAn et al. 1 120071 ) 



(120071) used 4 dwarf s to derive -1-0.11 ± 0.06, on a scale where the Hyades 



Pace et al. 



(120081) obt ain +0.27 + 0.10 from 7 dwarfs. As a final note 



More recently, lAn et al 
has [Fe/H] = -^0.13, while 
to the summary from the literature, I Taylor I ([200a) has attempted a rigorous assessment of the red- 
dening of the Hyades, Coma, and Praesepe and concludes that while the former two clusters have 
effectively no reddening, Praesepe has E(i? — V)= 0.027. His reevaluation of past high-dispersion 
work in this context implies [Fe/H] = -1-0.01 =b 0.04 for the cluster. 



With the exceptions of lTavlor I (|2006l ) and lPace et al. I (j2008l ) at opposite ends of the metallicity 
limits, there is general agreement that Praesepe is metal-rich compared to the Sun and similar to, 
or slightly lower than, the Hyades. This meshes well with the high-dispersion analysis of a few 
dozen Praesepe dwarfs by the Indiana group which, when processed on the same temperature scale 
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as NGC 3680 with E{B 
[Fe/H] = +0.15. 



V) = 0.0, leads to [Fe/H] = +0.144 on a scale where the Hyades has 



The results for the combined Li data for the Hyades and Praesepe (jSteinhauer 1 12003| ) are 
shown as a function of — ^ in Fig. 8. Blue circles and red squares are the Hyades and Praesepe 
data with measured abundances, respectively, while t he blue and red t r iangle s are upper limits 



The o r iginal EW data for Li i n the H yades are from iThorburn et al. I (|l993l ): iBoesgaard et al 



(ll988ll:lBoesgaard Trioicco I (Il986all: 



Soderblom et al. 



(Il990l). while the Praesep e data are from 



Boeseaard fc Budge I (|l988l ): ISoderblom et al. I (|l993l 'l: iBurkhart fc Coupnn ^1998 ). The EW data 
have been convert ed to A (Li ) using the same approach ad opted for NGC 3680. The B — V colors for 
Praesepe are from I Johnson I (jl952l ) and iMendoza I (|l967l ): both sets of photometry ar e in excellent 



agreem ent with each othe r. For the Hyades, we have adopted the photometry of iJoner et al. 



(1200a) whenever possible. Ijoner et al. I (120061) identify the need for a z ero-point shift of +0.008 
mag for the more commonly used data f r om I Johnson &: Knuckles I (j 19551 ): for stars for which only 
the photometry of I Johnson &: Knuckles I (jl955l ) is available, this offset has been applied. 



Though it is by no means perfect, the agreement between the superposed data is impressive, 
especially in the region of the Li-dip. The one obvious disagreement is in the B — V region between 
0.55 and 0.6. Praesepe exhibits a population of stars that has clearly depleted the surface Li by a 
significant amount compared to the well-defined upper bound in this color range. The depletion is 
non-uniform among the Praesepe stars, with one star reaching a factor of 20 below the expected 
value, and it is not present in the less populous sample from the Hyades. In the region of the Li-dip, 
however, we have drawn by eye a mean relation through the combined sample, as illustrated by the 
solid line. Within the center of the Li-dip, the points defining measured Li abundances have been 
given significantly more weight than the points with only upper limits. Three conclusions emerge 
from the mean relation: (a) the B — V center of the Li-dip is the same for both clusters at 0.44 ± 
0.01. For this metallicity and age, the isochrones predict that the stars at this point will have a 
mass of 1.421 ± 0.015 Mq] (b) the high mass limit of the dip occurs precipitously at B — V = 0.405 
± 0.005, equivalent to 1.494 ± 0.008 Mq while the low mass boundary occurs at B — V = 0.517 
and a mass of 1.28 Mq; and (c) the Li-dip is not symmetric. In particular, there is evidence that 
in the color range from B — V = 0.47 to 0.52, Li depletion has occurred at a higher rate than for 
stars just redward of this region, creating a secondary dip among these stars on the main sequence. 

If we assume that the Li-dip profile as defined by the Hyades and Praesepe represents that of 
IC 4651 at the same age, we can now translate this relation to the age, distance, and reddening of IC 
4651 and see how they they compare. The initial step is the evolution of the Li-dip to the age of the 
cluster. This is done by converting the Li-{B — V) relation to a Li-Mass relation, then converting 
it to a Li- My relation at the age of IC 4651. With the apparent modulus known, this readily 
translates into a predicted location for the Li-dip in the CMD turnoff region. The expectation is 
that the center of the dip should occur at the appropriate magnitude but, at best, the IC 4651 
data should lie at or below the mean relation at a given V due to potential depleti on in Li between 



an age of 0.7 Gyr and 1.5 Gyr. The Li data for IC 4651 come from three sources. iPasquini et al 
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(|200J), iBalachandran. Anthony- Twarog &: Twarog I (|l99ll ) and unpu blished data from Suc hitra 



Balachandran (private communication). Of these several sources, onlv iPasquini et al. I (|2004l ) lists 



measured equivalent widths as well as final abunda nces and stellar temperatures. We employed 



(B — V) colors from lAnthonv-Twarog et al. I (|l988l ) to compute temperatures for the dwarf stars 



consistent with our analysis path for NGC 3680 and NGC 7 52, using a redd e ning value E(i? 



V) = 0.12 and an assumed metallicity of [Fe/H] = 0.15. As iPasquini et al. I (|2004l ) list lithium 
abundances derived for two separate temperatures, one based on photometric colors, the other 
based on spectroscopic reductions, we were able to adjust the published Li abundances to conform 
with our temperature scale by applying a correction AA(Li) = 7.1 x 10~^AT. Results for stars 
common to more than one source were averaged if both results were derived abundances; if either 
or both results were upper limits, the more stringent upper limit was adopted. 

Initial attempts to superpose the predicted Li-dip distribution upon the observed trend in IC 
4651 made use of the Padova isochrones and failed in that the predicted location of the Li-dip at 
1.5 Gyr was 0.4 mag fainter than observed, implying that the mass center of the Li-dip for the 
Hyades and for IC 4651 were different. Part of the issue may be tied to the need to adjust the 
Padova isochrones to our adopted solar color and luminosity, shifting the color for a star on the 
main sequence at a given mass. While the combined shift was tested for the Hyades and confirmed 
in that it gave the appropriate distance modulus, it may well be that the relative required combined 
shifts in B — V and V are different at higher [Fe/H] from what they are at solar metallicity, thereby 
altering the mass-color relation. The switch to the isochrones was made, in part, because the 
offsets applied to the isochrones are smaller and conceivably would be less affected by alterations at 
metallicities different from solar. This has proven to be the case. Fig. 9 shows the Li distribution 
for the stars in IC 4651 (open circles are detections while triangles are upper limits) with the 
Hyades/Praesepe Li-dip relation evolved to an age of 1.5 Gyr and shifted by 10.25 mag, identifying 
the center of the Li-dip as V ^ 13.35. The mass of the stars at this location is 1.47 Mq, larger 
than the Hyades-based 1.42 Mq and a predicted Li-dip location V = 13.52. Note that a star 
with 1.47 Mq would have B — V = 0.416 at the metallicity and age of the Hyades, only slightly 
bluer than our adopted Li-dip center of B — V = 0.44. 

Given the small number of stars used to define the IC 4651 profile, the consistency between the 
shift predicted by the isochrones and the observations is encouraging. On the brighter, high-mass 
side of the dip, there appears to be no decline from A(Li) values near the solar system value until one 
reaches the stars evolving through the hydrogen-exhaustion phase and toward the red giant stage. 
Within the broad uncertainties, the sharp drop seen on the high-mass side remains unchanged. 
On the low-mass, fainter side, all the cluster stars fall on or below the Li-dip, as expected if some 
depletion has occurred from the upper limits for these stars as defined by the Hyades relation. Note 
also that the secondary dip identified within the Hyades will be difficult to verify due to the small 
sample populating the correct magnitude range. 
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5. Lithium in Open Clusters 

5.1. Li-Dip: Structure and Evolution 

To optimize the insight gained from the ehisters, we combine Fig. 7 and Fig. 9 in Fig. 10. 
Because the metallicities of NGC 3680 and NGC 752 are so similar and the age difference partially 
compensates for any evolutionary effect, no adjustment is made to the relative positions of the 
points (blue) in Fig. 7. The combined samples for these two clusters are at an effective true 
modulus of (m — M)^ = 10.12; if the metallicity is lowered to [Fe/H] = -0.12, the true modulus 
is reduced to 10.07. For IC 4651, the true modulus is (m — M)o = 10.03, identical within the 
uncertainties to the distance to NGC 3680. For consistency, however, we have added 0.09 mag to 
the V magnitudes for the data points in IC 4651 (red) and the mean relation (solid line) of Fig. 9 
to place them at the same distance as the other clusters in Fig. 10. 

Clearly, the Li-dip profile as defined by the combined Hyades and Praesepe data and evolved 
to an age of --^1.5 Gyr does an excellent job of mapping out the region of the observed composite 
Li-dip. Several key features are enhanced in the composite sample. (1) The upper bound at A(Li) 
= 3.35 for the stars above the Li-dip is well established. One star at the bright edge of the Li-dip 
sits above the upper bound by only 0.1 dex, an insignificant offset given the errors; (2) Keeping 
in mind that non-interacting binaries will be shifted horizontally to the right and that the sample 
of stars above the Li-dip is incomplete, the upper bound appears to end abruptly at V = 12.35. 
Stars brighter than this level arc leaving the turnoff for the giant branch, cither in or beyond the 
hydrogen-exhaustion phase. It should be emphasized, however, that not every star within the main 
sequence red hook exhibits reduced Li; (3) The vertical drop in A(Li) on the high mass side of the 
dip occurs over a very small range in V, consistent with the sharp Hyades/Praesepe drop with color. 
Likewise, the profile on the low-mass side matches the observational upper bound extremely well. 
This appears to indicate that over the ~1 Gyr of evolution between the Hyades and the average age 
of these clusters, little significant alteration has occurred in the Li distribution within the Li-dip, 
implying that the physical process producing the dip may have run its course by 0.7 Gyr. It should 
be remembered, however, that due to the steepness of the profile and the logarithmic basis of the 
plot, a drop in Li by a factor of 2 or less would be difficult to detect; (4) With the exception of a 
handful of stars near V = 14.7 in NGC 752, the average star fainter than the Li-dip is typically 
0.25 dex lower in A(Li) than the mean relation for the Hyades, though this differential would be 
cut in half if the mean dip relation were shifted to brighter V magnitudes by 0.15 mag. Within 
the errors, there does appear to be a modest decline in A(Li) for the cooler dwarfs outside the dip 
between the ages of the 0.7 and 1.5 Gyr. 

To close our discussion of the Li-dip, we can estimate the mass dependence of the center of the 
dip on metallicity. Our comparison of the Hyades/Praesepe dip (M = 1.42 Mq) with the evolved 
location in IC 4651 (M = 1.47 Mq) implies that, due to the uncertainties in the evolutionary 
models, independent of the uncertainty in the exact placement of the center of the dip for the older 
clusters, it is possible that the isochrones will predict too high a mass when matched to the CMD 
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position of the clusters near 1.5 Gyr. For NGC 3680 and NGC 752 at an approximate [Fe/H] = 
-0.07, the mass of the Li-dip center is 1.35 Mq. Contrasting this with IC 4651, a hnear relation 
over this metallicity range would be Mdip/MQ = 1.38(±0.04) + 0.4(ib0.2)[Fe/H]. For comparison 



using field stars over a much wider range in [Fe/H] and a different set of stellar models. IChen et al. 
(|200lh found Mdip/MQ = 1.42 + 0.28[Fe/H]. 



What is the physica l origin of the Li dip and the regions surrounding it? Once the Hyades 
and field dwarf Li data of Boesgaard &: Tripicco ( 1986al lbl) firmly established t hat F stars depleted 
far more Li than predicted by standard theory (see discussion in iDelivannis I (|200d )). a variety of 
mechanisms beyond the SSET were proposed. These loosely fall into three categories: mass loss, 
diffusion, and slow mixing due to rotation or wa ves. In the mass loss scenario, the sm all outer 
mass fraction where Li is preserved is simpl y lost (ISchramm. Steigman &: Dearborn 1119901 ). In the 



original diffusion scenario for the Hyades (|Michaud 1 119^ diffusion is unimportant for surface 



Li depletion in G dwarfs, but as one goes to the hotter early F dwarfs, the increasingly shallow 
SCZ allows increasingly efficient draining of the Li out of the base of the SCZ via gravitational 
settling and thermal diffusion. This creates the red side of the gap. In the middle of the Li gap, Li 
retains an electron and with its now much larger cross-section, Li becomes radiatively accelerated 
upward, creating the hot side of the gap. However, this can cause large Li overabundances in 
stars hotter than B — V =0.40 (a "Li peak"), which was not observed, so a finely-tuned mass- loss 
rate was included i r i such stars to balance the diffusionj^. The more refined diffusion models of 
Richer &: Michaud I (|l993l ) predicted a much narrower Li dip (~100K in Te), which could no longer 
explain the full width of the observed Li dip (~ 300K for the primary dip and 500K or more when 
including the secondary dip), and even less so the more widespread Li depletion hotter and cooler 
than the dip. 



Mechanisms that can drive sl ow mixing include, a) waves caused by the SCZ (jGarcia Lopez &; Spruit 



1991 



Talon Charbonnel II2003I ). and b) rotation-related effects such as m eridional circulation or 



the on s et of angular momentum transport induced primarily by secu lar shear (jPinsonneault. Kawaler Demarque 



(|l990l ): IChaboyer et al.l (1199511: IDelivannis &: Pinsonneault I (|l997l ). collectively referred to as the 
Yale models). Charbonneau &: Michaud ( 19881 ) proposed that meridional circulation combined 

with diffusion could explain the Li dip, though later work suggested that the efficiencv of meridional 

circu lation must be greatly suppressed (jCharbonneau. Michaud &: Proffitt 1ll989l :l iPinsonneault. Kawaler &: Demarc 
1990) • Another candidate process, turbul ent transport, overdepletes stars beyond the hot side of 



the Li dip when diffusion is allowed for (jCharbonneau &: Michaud Ill988l ). In the Yale models. 



surface angular momentum loss acting primarily during the early main sequence of G dwarfs exac- 
erbates differential rotation with depth, which sets off a secular shear instability, causing mixing. 
Different initial rotation rates cause the observed dispersions in Li abundances. Increasing initial 



^Li peak stars may exist afterall: IDelivannis. Steinhauer fc Jeffries I ( 2003 ) discovered a star in tlie Hyades-aged 
cluster NGC 6633 near 7100K with A(Li) = 4.3, at least an or der of magnitude above the assumed initial abundance of 
the cluster. This i s consistent with the diffusion predictions of Richer fc Michaud 1 1993h . although other explanations 
are also possible ( Laws fc Gonzalez 20031 : Ashwell et al. 20051 ) 
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angular momentum with mass and corresponding greater mixing results in the cool side of the dip. 
But stars progressively hotter than the break in the Kraft curve (which coincides with the middle 
of the Li dip) lose progressively less angular momentum, thus producing the hot side of the Li dip. 

These various classes of models meet with varying degrees of success in reproducing the Hyades 
Li dip, and distinguishing between them using the Hyades dip alone can be challenging. Fortunately, 
there exist additional powerful discriminators, that include, a) the timing of the formation of the Li 
dip and the evolution of its Li-mass morphology, b) the Li-Be and Be-B depletion correlations, c) 
short-period-tidally-locked binaries, and d) subgiants evolving out of the Li dip, whose deepening 
convection zones reveal the shape and depth of the Li preservation region of their main sequence 
progenitors. 

The early onset of Li depletion by factors of 4 to 6 and over a wide T^ range (6000 to 6700K, 
at minimum) observed in the 150-Myr-old cluster M35 (jSteinhauer &: Delivannis 1 12004 ) favors the 
Yale rotational models as the dominant effect at early ages and argues against mass loss or diffusion 
which require longer time scal es. The striking Li-Be depletion correlation disc overed in field dwarfs 



on the cool side of the L i dip (|Delivannis et al. Ill998l : iBoesgaard et al. Il200ll ) and in open clusters 



(|Boesgaard et al. II2004I ) argues strongly in favor of slow mixing induce d by rotation, as in t he Ya le 
and other models that produce the specifically observed Li/Be ratio (e.g lCharbonnel et al. I ()1994l )). 
and argues strongly against diffusion (where Li and Be are depleted from the surface at a similar 
rate) and mass loss (where all the Li must disappear before Be depletion can be observed). The Li- 
Be correlatio n can even discriminate among models with slow mixing; for example, the wave-driven 
mixing in the lGarcia Lopez &: Spruit I (jl99ll ) models does not penet rate very deeply, and a s a res ult 
the predicted Li/Be ratio is too large. The wave-driven mixing of iTalon &: Charbonnel I (|2003l ) is 
more promising, especially in its association of angular momentum transport by waves, and in its 
potential to explain the spread of Li abundances as a function of Tg for stars cooler than the Li dip, 
but it remains to be seen whether more quantitative models can reproduce constraints other than 
those from Li in the Hyades alone (such as, but not limited to, the constraints described here). Since 
Be is preserved to roughly twice the depth as Li with B preserved to a depth twice as large again, the 
observed Be-B depletion correlation furthe r illustrates the depth t o which mixing occurs, and further 
defines its r elative efficiencv w i th de pth (IChaboyer et al.lll995l ). Combining the Yale rotational 
models with IZahn &: Bouchet I (|l989l )'s tidal circularization theory results in the prediction that 
tidally locked binaries with sufficiently short period should exhibit higher Li abundances than 
single stars because they exchanged angular momentum during the early pre-main sequence when 
their interiors were not yet hot enough to destory Li. Thus, they do not exhibit the subsequent 
mixing and Li dep l etion induced by angular momentum loss that normal, single stars do. (See also 



Soderblom et al. I (|l990l ).l This prediction is supported by the observatio n of high-Li short perio d 



binaries above the late-F /early G Li plateau and elsewhere in the H yades (IThorburn et al. Ill993l ) 



in M67 (jPeliyannis et al. Ill994l ). and in a variety of other contexts (iRyan &: Deliyannis Ill995l ). Li 



in subgiants in the solar-aged M67 cluster reveal that the Li preservation region during those stars' 
MS evolution was smaller that that predicted by standard or diffusion models, consistent with 
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that predicted by the Yale rotational models, and larger than that predicted by mass loss. These 
considerably varied diagnostics all consistently favor slow mixing (probably induced by rotation) as 
the primary cause of Li depletion in F dwarfs, and consistently argue against mass loss or diffusion. 

The present study adds additional interesting constraints that successful mechanisms must 
meet. The possible delineation of the Li dip into a deep primary and shallower secondary Li dip in 
the Hyades+Praesepe data, which is also consistent with the Li data in the older clusters NGC3680, 
NGC 752, and IC4651, suggests that Li depletion occurs fairly early (before 200Myr), over a wide 
range of Tg (at least 6000-6700K) and is caused primarily by rotationally-induced mixing, but a 
more severe Li depletion continues to develop in a more restricted Tg range near 6700K later but by 
the age of the Hyades. This further depletion could be caused by or be aided by diffusion. Although 
two stars with very low Li in the p rimary Hyades Li dip are also known to have a significant Be 



depletion (jBoesgaard King II2002I ). consistent with rotational mixing, the Li data are only upper 
limits, and thus the Li/Be is not known in these stars, and could allow for the possibility of late Li 
diffusion, especially at a Te where the picture gets furt her complicated by the possibility of upward 



radiative acceleration of Be (jRicher &: Michaud Ill993l ). even as Li continues to diffuse downward. 



Additional Li data, preferably detections, in the deep parts of the primary Li dip, ideally along 
with Be data, could help cast light on whether diffusion ever contributes to the formation of the 
Li dip. In NGC 3680 stars hotter than the Li dip, we stress that we have found a Li spread of at 
least 1 dex in magnitude. That is, these stars do not exhibit a constant, presumed undisturbed 
initial cluster Li abundance. (However, a larger fraction of IC 4651 stars, indeed most of them, are 
found near the Li upper bound; only three stars have been observed in NGC 752, but they also lie 
close to the upper bound.) Such a spread in NGC 3680 is not altogether surprising, as the SCZs 
are extremely shallow in these stars, and many complex competing factors can affect the surface Li 
abundances, including rotationally-induce mixing of varying types, diffusion, magnetic fields, and 
even just tiny rates of mass loss. 



5.2. Li in the Giants 



While the focus of our study has been the turnoff dwarfs, the ad ded data on the giants in 
NGC 3680, as well as additional data in IC 4651 (jPasquini et al. II2004I). promise poten tial insight 
into post-main-sequence evolution of Li as discussed by ATT and lPasquini et al. I (j2004 ). For NGC 
3680, the sample of PRP has been increased from six to nine giants; all three members of the 
expanded sample are probable binaries. As emphasized in ATT, the surprising result from the 
earlier work of PRP was the discovery that while the Li abundances among the giants of NGC 3680 
were reduced relative to the turnoff stars, of the six giants studied, four had measurable abundances 
rather than upper limits, and three of the four measurable values sat near A(Li) = 1.0. This led 
to a bimodal distribution in A(Li). The combined sample confirms and expands the result. Of the 
three new giants, 2 have A(Li) = 0.98 and 1.09, while the third has 0.50. For the 9 giants, five have 
an average A(Li) of 1.08 ± 0.08, while the remaining 4 have A(Li) = 0.45 it 0.07 or an upper limit 
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measurement. It should be noted that the true hmit for KGP 1379 is undoubtedly lower since the 
unusual blue color of this star implies a likely main sequence turnoff companion. While the Li line 
strength would be weakly impacted by the companion, the effect of the bluer color is to raise the 
effective temperature, producing a significant overestimate for Li. 

How does NGC 3680 compare to NGC 752? First, of the three clusters, NGC 3680 is clearly the 
oldest, while NGC 752 and IC 4651 are virtually identical in age if one accepts the higher reddening 
value for the latter cluster. As discussed in ATT, NGC 752 and NGC 3680 have distinctly different 
A(Li) distributions among the giants. For NGC 752, of 11 highly probable members on the giant 
branch, only 2 have detectable Li abundances, both above A(Li) = 1. The remaining nine stars 
have only upper limits of 0.5 or less. 

Beyond the age difference, the second way in which NGC 752 and NGC 3680 distinctly differ 
is in the CMD distribution of the giants. A comparison of Fig. 3 with Fig. 5 shows that the 
giants in NGC 3680 sit squarely on the first-ascent giant branch in the luminosity range associated 
with the red giant bump. By contrast, while the giants in NGC 752 occupy a comparable range in 
luminosity, they sit systematically blueward of the first-ascent giant branch and 0.07 mag bluer in 
B — V than the typical star in NGC 3680, occupying the CMD area commonly associated with the 
red giant clump and He-core-burning stars. Finally, the two stars in NGC 752 with measurable Li 
do not sit among the stars populating the clump, but lie along the red edge of the CMD distribution 
(see Fig. 8 of ATT). A qualitative explanation comes from the observation that some stars evolving 
through the red hook at the turnoff exhibit reduced Li from the limiting value of 3.35. Note that 
of all the stars with measured Li in the three clusters, only one star in IC 4651 is situated in the 
subgiant region of the CMD between the turnoff and the vertical first-ascent giant branch. It has 
A(Li) = 1.5. All other single stars with reduced Li populate the red ho ok. The implication is 



that some form of unexpected mixing, possibly due to enhanced rotation (jPasquini et al. II2004I ). 
is destroying the atmospheric Li as the star evolves away from the main sequence and up the giant 
branch. Thus, first-ascent red giants may still retain a detectable level of Li near A(Li) = 1 as they 
pass through the red giant bump phase and/or bring new Li to the surface as a byproduct of the 
interior evolution and mixing that creates the red giant bump. Evolution to the He-core-burning 
phase reduces A(Li) of 0.5 or lower, unless an alternative source of Li production is found. When 
coupled with the timescales in the different evolutionary phases, giant branch abundances should 
show a natural dichotomy. One group of predominantly first-ascent giants will have measurable Li 
near A(Li) = 1, while the more evolved stars will have A(Li) = 0.5 or less. NGC 752 is dominated 
by the latter group, NGC 3680 by the former. 

Given the pattern above, it is straightforward to predict where IC 4651 should fall in the Li- 
dichotomy. Its age is effectively identical to NGC 752 but the mass of a typical star on the giant 
branch (1.93 Mq) is higher than that for NGC 752 (1.86 Mq) or NGC 3680 (1.75 Mq). The stars 
on the giant branch, especially with the increased cluster reddening, sit squarely blueward of the 
first-ascent giant branch by more than 0.1 mag va. B — V (Fig. 7). The Li distribution among the 
giants should resemble that of NGC 752. As discussed in ATT using the unpublished data supplied 
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by Suchitra Balachandran, of 11 member giants observed in IC 4651, only one had detectable Li 
with A(Li) = 1. The other 10 star s only had upper hmits, with the limits ranging from A(Li) = 



1.05 to -0.1. iPasquini et al. I (|2004l ) observed 5 giants, three of which overlapped with the earher 
sample. For the two new giants, the stars had upper limits to A(Li) near 0.2. For the stars in 
common to the surveys, the newer data showed an upper limit of 0.38 where the older sample found 
an upper limit of 0.2, a measured value of 0.5 where the older data implied an upper limit of 1.0, 
and a measured value of 0.82 where the previous result was an upper limit of 0.75. Therefore, of 
13 observed stars, one giant definitely has A(Li) at 1.0 with possibly a second near 0.8. Eight more 
giants have a measured value of 0.5 or an upper limit of 0.7 or less. The remaining three giants 
have upper limits, but they range between 0.95 and 1.05, so the possibility of detectable Li cannot 
be ruled out. If we class these last three as unknowns, the fraction of giants with high Li (A(Li) 



1) in IC 4651 is 2 out of 10, statistically the same as in NGC 752. 



6. SUMMARY 

The specific focus of this study has been the Li abundance distribution of stars at the turnoff of 
the intermediate-age open cluster, NGC 3680, with the general goals of linking future Li analyses 
of turnoff stars in NGC 6253 to the larger array of published data in other open clusters while 
reevaluating NGC 3680 itself in the context of the evolution of the Li-dip. After taking into account 
the modest differences in the reddening and temperature scales, there is very good agreement with 
past Li abundance determinations. The combined data set more than doubles the turnoff sample 
and increases the red giant sample from 6 to 9 stars. The expanded sample better delineates the 
Li-dip region and reinforces the observed dominance of the giant branch by stars with A (Li) near 
1. 

Spectroscopic analysis of the red giants in NGC 3680 produces [Fe/H] = -0.17 ± 0.07 from 8 
stars. Previous giant star spectroscopy gives the same result if the spectra are analyzed with the 
same stellar parameters. The absolute scale can, however, be set at any value between [Fe/H] = 
0.00 and -0.30 by adopting microturbulence velocities between 1.3 km/sec and 1.9 km/sec. For the 
dwarfs, the scatter is a strong function of the S/N; from 5 single stars at the turnoff with S/N above 
110, [Fe/H] = -0.07 it 0.02. There is no statistically significant evidence from lines other than Fe 
that NGC 3680 has anything but scaled-solar elemental abundances. An average of all abundance 
determinations weighted by the inverse of the standard errors of the mean heavily favors the dwarfs 
and produces [Fe/H] = -0.08 it 0.02 (sem). With the reddening and metallicity known, isochrone 
comparisons require an age of 1.75 Gyr and an apparent modulus of (m — M) = 10.3. These data 
place the center of the Li-dip at a 1.35 Mq. 

With the fundamental cluster information in hand, NGC 3680 is contrasted with two slightly 
younger clusters, NGC 752 and IC 4651, and the even younger clusters, Hyades and Praesepe. 
Assuming a slightly higher reddening value of E(i? — V) = 0.12 for IC 4651, with [Fe/H] = - 
0.06 and -1-0.15 for NGC 752 and IC 4651, respectively, both clusters have an age near 1.5 Gyr, 



-32- 



with apparent moduli of 8.42 and 10.4, respectively. For the Li-dip comparisons, the Hyades and 
Praesepe are both assumed to have zero reddening and [Fe/H] = +0.15; the distance moduli and 
any small age difference between these two clusters are irrelevant for purposes of the discussion. 
The combined Hyades/Praesepe samples are used to delineate a Li-dip profile, centered at B — V 
= 0.44 or 1.42 Mq. Evolution of the profile to the mean age of the three older clusters generates a 
CMD morphology that reproduces the observed Li-dip turnoff extremely well, implying that while 
the mass of stars in the center of the dip varies positively with metallicity, the color profile of 
Li-depletion across the dip may be similar but with the central color shifted to bluer values for 
lower [Fe/H]. Over the range from [Fc/H] = -0.1 to +0.15, the central mass of the dip follows the 
relation M/M0 = 1.38 (±0.04) + 0.4 (±0.2)[Fe/H]. Within the errors, there appears to be little 
change in the degree of Li-depletion for stars within the dip for ages greater than the Hyades. For 
lower mass stars just below the Li-dip, there is some evidence for modest depletion in Li (^0.2 dex) 
between 0.7 Gyr and 1.5 Gyr, assuming the clusters all formed with the same level of Li on the 
ZAMS. For the giants, the color and Li dichotomy discussed in previous investigations has been 
confirmed and strengthened for NGC 752 and NGC 3680. The younger cluster has, relative to the 
first-ascent giant branch, a significantly bluer giant branch distribution than NGC 3680, while only 
2 of 11 stars have detectable Li falling within the A(Li) = 1 bin, in contrast with 5 of 9 for NGC 
3680. For IC 4651, the comparable numbers are 2 of 10, as expected given its similarity to NGC 
752, but these claims are tempered by the uncertainty in the true Li abundance for 3 additional 
giants. If this pattern holds up, it strengthens the argument that, due to the lower mass of the stars 
leaving the main sequence, the giants populating the supposed clump in NGC 3680 are first-ascent 
giants, placing an important constraint on the transition in evolved structure for stars in the mass 
range of 1.8 to 1.6 Mq. 

The three clusters profiled in this investigation are regularly compared because they are so 
similar but also because they represent objects that, by the standards of open cluster work, are well- 
studied. The spectroscopic sample of stars within NGC 3680 has been doubled by this analysis but 
the Li-dip profile, as presented in Fig. 10, can only be adequately mapped through the superposition 
of data from all three clusters. One reason for the frequently qualified conclusions in the present 
paper, as well as others related to Li evolution, is the small number of single members, as is the 
case in NGC 3680; in other clusters, the paucity of data represents an inadequate fraction of cluster 
members incorporated in analyses. As an example, the plateau effect for stars above the Li-dip 
seems apparent in Fig. 10, but delineation of the brighter stars' depletion as they evolve toward 
the giant branch is hampered by observations of only one star in this region in NGC 752 and only 
five of more than two dozen possible candidates in IC 4651. A more complete set of observations 
has been obtained in NGC 752 by one of the authors (CD), but an expanded sample in IC 4651, 
especially among the giants, could prove invaluable. 

In the arena of future observations, two needs come to mind. One of the weaknesses of 
the mctallicity-dcpendent relation derived for the mass of the center of the Li-dip is the small 
baseline in [Fe/H]. A critical test of its reliability could come from the observation of a more metal- 
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poor open cluster in a comparable age range. The obvious options are NGC 2420 and NGC 2506 



2006 



with typical estimated abundances (jTwarog et al. 1119971 : iFriel et al. II2002I : lAnthonv-Twarog et al 



Lee et al. 1120081) b e tween [Fe/H] = —0.3 and —0.4, though, for a higher estimate for NGC 



2506, see lCarretta et al. I (|2004l ). The second somewhat surprising deficiency is the lack of modern, 
broad-band, wide-field CCD photom etry for two of t he three clusters discussed, IC 4651 and NGC 
752. The uvby survey of IC 4651 by iMeibom I (l200d ) is outstanding, but a complementary broad- 
band data base would make isochrone comparisons more straightforward, particularly on the lower 
main sequence; the composite sample for NGC 752 clearly has a high but uneven degree of precision 
and remains incomplete at fainter magnitudes. 
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Fig. 1. — Proper- motion probability distribution for stars identified as most likely to be potential 
members of NGC 3680. 

Fig. 2. — Radial velocities for the stars in the current study as a function of the proper- motion 
probability. The dashed horizontal line defines the cluster mean velocity from the probable dwarf 
members, while the horizontal dash-dot lines define the radial-velocity limits for probable single-star 
cluster members, between —1.8 and 5.1 km/sec. 

Fig. 3. — CMD for stars observed in the present study, excluding probable non-members. Filled 

circles are likely cluster members while open circles identify members that are suspected of being 
binaries. Triangles denote stars observed by PRP, with binaries similarly denoted by open symbols. 
Superposed are isochrones of appropriate metallicity with ages of 1.7 Gyr and 1.8 Gyr, adjusted 
for the reddening and apparent distance modulus of NGC 3680. 

Fig. 4. — A(Li) vs. V for stars in NGC 3680. Filled circles and triangles are stars observed in 
this survey with detectable Li and upper limits, respectively, combined with the data of PRP if 
the samples overlap. Open circles and triangles are stars with detectable Li and upper limits, 
respectively, observed only by PRP but transferred to our temperature and Li scale. 

Fig. 5. — CMD for NGC 752 superposed upon a isochrones with [Fe/H] = -0.05, (m - M) = 
8.4, and E{B — V) = 0.03 and ages of 1.4 Gyr and 1.5 Gyr. Known binaries and non-members have 
been removed from the CMD. 

Fig. 6. — The combined Li data for NGC 3680 (black symbols) and NGC 752 (red symbols). Circles 
are detected values while triangles are upper limits. 

Fig. 7. — CMD for IC 4651 superposed upon a 1.5 Gyr isochrone with [Fe/H] = +0.15, (m — M) 
= 10.4, and E{B — V) = 0.12. Triangles are stars identified as spectroscopic binaries. Known 
non- members have been removed from the CMD. 

Fig. 8. — A(Li) as a function of B — V for stars in the Hyadcs (blue symbols) and Praescpc (red 
symbols). Circles are detectable values of Li while the triangles are upper limits. Solid line is an 
estimate of the mean relation through the points in the region of the Li-dip. 

Fig. 9. — A(Li) as a function of V for stars in IC 4651. Circles are detectable measures while 
triangles are upper limits. The solid line is the mean relation of Fig. 8 evolved to the age of IC 
4651 and shifted in V by 10.25 mag. 

Fig. 10. — A(Li) as a function of V for the combined cluster sample of Fig. 6 and Fig. 9. Blue 
circles are detectable measures while triangles are upper limits for NGC 3680 and NGC 752; the 
red symbols have the same meaning for IC 4651. The solid line is the mean relation of Fig. 9. The 
solid relation and the data for IC 4651 have been shifted to fainter magnitudes by 0.09 mag in V. 
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Table 1. Basic Data for Program Stars in NGC 3680 
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0.6 


PhotBin(M95;AHTC) 


1395 


37 


37 


13.030 


0.433 


53 


22.2 


0.7 


1.8 


1.0 




1404 


27 


27 


10.770 


1.134 


71 


16.0 


0.3 


2.2 


0.6 


SBlO(M95,N97,M07) 


1405 


31 


31 


13.135 


0.456 


69 


21.6 


0.6 


0.4 


0.8 


SB1(N97) 


1410A 


240B 




12.424 


0.437 




60.3 


3.1 


1.1 


4.1 


SB?(KP97); w/sb? 


1446 


2058 




15.801 


0.795 


16 


10.6 


0.4 


51.2 


0.9 




1454 


32 


32 


12.917 


0.478 


66 


21.2 


0.5 


2.0 


0.8 


SB1(N97) 


1461 


41 


41 


10.930 


1.121 


23 


16.6 


0.3 


1.0 


0.5 




1469 


44 


44 


10.002 


1.235 


67 


17.2 


0.3 


1.2 


0.6 


SB? (KP97) 


1506 


42 


42 


12.977 


0.441 


62 


51.4 


2.8 


-2.7 


3.7 


sb? 


1507 


46 


46 


13.358 


0.452 


60 


16.0 


0.4 


2.0 


0.7 




1510 


43 


43 


12.552 


0.448 


68 


30.6 


1.6 


2.6 


1.9 


PhotBin(ATTS) 


1522 


45 


45 


12.034 


0.504 


67 


30.6 


1.1 


1.2 


1.4 




1524 


82 


46A 


13.498 


0.485 


49 


29.3 


0.8 


1.4 


1.0 


SB?(KP97) 


1534 


1085 




14.309 


0.565 


42 


5.6 


0.1 


2.6 


0.6 




1554 


1124 




14.285 
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51 
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15.404 


0.712 


13 
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0.3 


-0.1 
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1624 


2062 




13.634 


0.467 


61 


18.3 


0.4 


-1.8 


0.7 


SB1(N97) 
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0.867 
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51 
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0.5 


4.4 
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1019 
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1815 


1076 




15.241 


0.711 


6 


11.2 


0.3 


-12.9 


0.6 




1873 


53 


53 


10.859 


1.120 


69 


18.0 


0.3 


2.0 


0.7 





Table 2. Lithium EW and Abundance Data for Program Stars in NGC 3680 



KGP Temp. EW EW' gew A(Li) ua S/N per pixel FWHM 



Abundances from detected Lithium Lines 
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Table 2 — Continued 



KGP Temp. EW EW gew A(Li) a a S/N per pixel FWHM 



Upper Limit Lithium Abundances 
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Table 3. Metal Abundance Data for Program Stars in NGC 3680 



KGP 


[Fc/II] 






[Si/II] 


[CVi/II] 


[Xi/II] 






Giants 


1175 


-0.13 


0.13 


8 


0.20 




-0.15 


0.04 


1261 


-0.18 


0.11 


8 


0.00 




-0.33 


0.05 


1286 


-0.22 


0.06 


8 


-0.16 
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0.07 
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0.12 
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0.00 




-0.21 


0.20 


1404 


-0.10 


0.13 
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0.04 
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0.09 


8 


0.09 
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0.06 
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-0.14 


0.16 


8 


0.17 




-0.09 


0.04 


1873 


-0.09 


0.23 


8 


-0.02 
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0.11 
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0.00 
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6 


0.02 


-0.05 


-0.08 


0.26 
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0.22 


6 
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0.37 
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0.17 


6 
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0.15 
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0.05 


0.17 


5 


0.04 
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0.00 


0.32 
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0.15 


6 


0.04 
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0.07 


4 
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Table 4. Abundance Data by line for Program Stars in NGC 3680 



A [A/H] a [A/H] a N^t^rs 
Giants Dwarfs 



Iron 



6581.2 


-0.23 


0.16 
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-0.25 


0.11 
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0.09 


11 
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11 


Calcium 


6717.7 






-0.17 


0.12 


13 
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0.11 
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13 
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-0.20 


0.11 


0.03 


0.14 


13 
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0.11 
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0.16 


13 



Silicon 



6721.8 0.05 0.12 -0.06 0.10 



13 



